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FOREWORD 

This  Technical  Report  was  prepared  by  the  Power  Systems  Branch  of 
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Southall  and  Frederick  C.  Brockhurst,  Project  Engineers.  This  work  took 
place  over  the  period  from  1  Feb  1978  through  31  March  1979. 
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SECTION  I 
INTRODUCTION 

1.1  BACKGROUND 

Prediction  of  weights,  volumes,  and  efficiencies  for  components  in 
power  supplies  for  high-power  (multi -megawatt)  applications  is  an 
important  first  step  in  determining  feasibility  for  use  aboard  aircraft 
with  strict  size  and  weight  constraints.  Such  a  power  supply  system 
is  illustrated  in  Figure  1. 

In  1975  and  1976,  the  U.S.  Air  Force  Aero  Propulsion  Laboratory 
(AFAPL)  sponsored  several  contractual  programs  under  the  name  of  the 
“High-Power  Study"  which  resulted  in  computer  programs  predicting  size, 
weight,  and  performance  for  every  component  of  a  power  supply  system 
like  the  one  shown  in  Figure  1.  Programs  were  developed  for  turbines 
(References  1,2),  electrical  generators  (including  both  conventional, 
iron-core,  round  rotor  (Reference  2),  and  superconducting  (References  2,  3) 
and  power  conditioning  components  (Reference  4).  Some  of  these  programs 
were  detailed  design  programs,  while  others  used  curve  fits  to  detailed 
point  designs  within  the  scope  of  the  study  and  were  therefore  not  useful 
for  power  or  voltage  requirements  outside  the  valid  parameter  space  of 
the  study. 


POWER  POWER  POWER 


Figure  1.  Generalized  High-Power  Airborne  System 
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Since  these  studies  were  completed,  the  AFAPL  has  begun  two  develop¬ 
ment  programs  in  advanced  technology,  second  generation  airborne  electrical 
generators  which  are  much  different  from  the  machines  considered  in  the 
High-Power  Study.  One  program  is  for  the  development  of  a  prototype  5  MW 
permanent  magnet  generator  (PMG)  which  would  use  high-energy-density 
samarium  cobalt  magnets  for  the  field  (Reference  5).  The  other  program 
is  a  20  MW  superconducting  generator  using  superconducting  wire  for  the 
field  coils  and  a  conductive  environmental  shield  (References  6,  7,  8). 
Algorithms  developed  in  this  report  are  specifically  related  to  these  two 
classes  of  electrical  generators. 

1.2  PURPOSE  OF  REPORT 

This  report  presents  simplified  weight  and  volume  algorithms  for  the 
two  classes  of  advanced  electrical  generators  discussed  above  and  indicates 
how  system  weight  is  affected  by  the  type  of  generator  selected.  The 
algorithms  are  simplified  in  the  sense  that  they  are  not  detailed  design 
programs  which  consider  every  electrical,  thermal,  and  mechanical  aspect 
of  electromechanical  machine  design.  They  are  simplified,  "approximate 
design,"  programs  which  use  fundamental  sizing  relationships  of  electro¬ 
mechanical  machine  design  coupled  with  current  technology  limits  of  the 
design  parameters  pertinent  to  advanced  permanent  magnet  and  super¬ 
conducting  generators.  Estimates  of  weights  and  volumes  predicted  by  the 
simplified  algorithms  are  shown  to  be  within  10%  of  current  design 
estimates  of  generator  weights  and  volumes.  Of  great  importance  is  the 
fact  that  the  algorithms  can  be  used  for  any  values  of  power,  voltage, 
rotational  speed,  etc.,  which  makes  the  algorithms  very  flexible. 

In  addition,  this  report  describes  the  important  considerations  of 
system  efficiency,  component  efficiency,  and  the  propagation  of  efficiencies 
in  a  power  supply  system. 

Finally,  generator  and  power  supply  component  algorithms  are 
exercised  to  illustrate  how  the  results  are  used  in  power  supply  system 
design.  This  study  compares  weights  and  volumes  of  the  two  types  of 
generators  at  output  power  levels  of  approximately  10  and  20  megawatts. 

In  addition,  overall  system  weight  is  calculated  for  seven  different  point 
designs . 
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1.3  FUTURE  WEIGHT  AND  VOLUME  COMPUTER  PROGRAMS 

The  ultimate  goal  is  to  have  algorithms  which  give  accurate  weight 
and  volume  estimates  for  every  component  of  the  system  shown  in  Figure  1 
for  a  self-consistent  total  system  size  and  weight.  Indeed,  as  dis¬ 
cussed  in  Section  V,  this  is  mandatory  for  a  valid  comparison  of  alter¬ 
native  power  sources.  Self-consistency  means  that  all  interface  parameters 
between  components  must  be  taken  into  account.  For  example,  the  generator 
output  voltage  and  frequency  influence  the  size  of  the  power  conditioning 
subsystem.  Although  this  can  be  done  manually,  as  in  Section  V,  an 
interactive,  computer-aided  power  system  design  program  has  been  developed 
under  contract  for  this  purpose  (Reference  9).  This  program  allows  the 
interactive  execution  of  all  the  computer  programs  developed  under  the 
High-Power  Study  to  arrive  at  a  total  system  weight  and  volume.  Also, 
the  program  can  be  modified  to  accept  new  programs  for  advanced  components, 
such  as  the  permanent  magnet  generator,  or  to  accept  improved,  revised, 
or  new  algorithms  for  other  components. 
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SECTION  II 

POWER  SUPPLY  EFFICIENCY  CONSIDERATIONS 

Because  of  the  influence  on  required  prime  mover  and  electrical 
power,  it  is  important  to  analytically  determine  the  effects  of  com¬ 
ponent  efficiencies  on  the  power  supply  system.  Figure  2  illustrates 
a  power  supply  system  which  supplies  primary  loads  P^  and  PL,,  and 
auxiliary  electrical  power  PAUX’  It  should  be  noted  that  the  generator 
exciter  power  (Figure  1)  is  assumed  to  be  zero,  which  is  valid  for 
permanent  magnet  and  very  nearly  correct  for  superconducting  generators. 
The  two  primary  loads  are,  in  general,  rated  at  different  power  and 
voltage  levels.  If  the  number  two  transformer  were  used  to  supply  power 
to  the  number  two  primary  load,  the  system  would  appear  as  in  Figure  3. 
Component  efficiencies  are  indicated  above  the  block  representing  each 
component. 

Efficiency  is  defined  as 


n  Pout/,Pin’ 


(1) 


where  PQut  is  the  power  out  of  a  component  and  P . n  is  the  power  into  a 
component.  Heat  generated  within  each  component  is  given  by 

"loss  ■  pi n'pout  ■  ^  Pour  <““«> 


From  the  definition  of  efficiency  and  the  specification  of  the  required 
output  powers,  P^  and  P^.  the  auxiliary  electrical  power  requirement, 

,  the  power 


P^ux,  and  the  auxiliary  mechanical  power  requirement,  P^^ 


AUX 


requirement  at  any  intermediate  stage  within  the  system  can  be  determined 
if  all  the  component  efficiencies  are  known.  For  example,  P^  ,  the  input 
to  the  number  one  power  conditioner  can  be  found  by  dividing  P^  by  n-j ,  i.i 


P1  =  PLl/nl 


(3) 
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Powers  at  intermediate  stages  upstream  can  be  calculated  in  a  similar 
manner.  Since  the  required  output  power  of  the  generator,  PQ,  is  of 
primary  interest  in  this  report,  expressions  for  PG  are  given  for  the 
systems  shown  in  Figures  2  and  3.  For  the  system  with  transformer  T1 
supplying  both  primary  loads  (Figure  2). 


p  _  n2  PL1  *  n1  PL2  + 
G  "  nT1  n-,  n2 


,  (watts) 


and  for  the  system  where  transformer  T1  supplies  one  primary  load  and 
transformer  T2  supplies  the  other  primary  load  (Figure  3). 


The  mechanical  shaft  power  which  must  be  supplied  by  the  turbine,  for  the 
system  of  Figure  2,  is  given  by 


PTURB 


Pl1  * P|-2  ♦  e-m -  ♦  V"  aux  _  (MMS) 

T1GB1  nTl  nl  n2  nGB1  nG  nT2  nGB2 


(6) 


and  for  the  system  of  Figure  3  by 


TURB  nGBi  nG  nT1  n1  nGB^  nG  nT2  n2 


MECH  AUX 


(watts) 


Generator  efficiency,  nG>  Is  calculated  as  part  of  the  generator 
algorithms.  pjG^B  is  used  in  estimating  the  amount  of  prime  mover  fuel 
required,  since  fuel  consumption  is  directly  proportional  to  shaft  power. 
Fuel  weight  is  a  major  part  of  the  overall  system  weight. 
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SECTION  III 

PERMANENT  MAGNET  HIGH-POWER  GENERATORS 


3.1  BACKGROUND 

This  analysis  provides  weight  and  volume  estimates  for  iron-core 
alternators  with  permanent  magnet  field  poles.  The  term  "iron-core"  is 
used  to  distinguish  this  class  of  machines,  which  relies  upon  a  relatively 
small  air  gap  and  a  high  permeability  iron  path  for  the  circulation  of 
magnetic  flux  from  rotor  to  stator,  from  machines  which  utilize  relatively 
large  air  gaps  (air-core  machines).  The  field  structure  in  the  latter 
case  must  provide  tremendous  magnetomotive  forces  (mmf)  to  establish  the 
required  magnetic  flux  densities  in  the  stator.  Large  air  gap  machines 
generally  have  superconducti ng  field  windings  and  are  described  in 
Section  IV. 

3.2  WEIGHT  AND  VOLUME  EQUATIONS 

The  primary  components  of  the  alternator  are  indicated  in  Figure  4. 
They  include  the  yoke  (or  back-iron),  the  teeth,  the  copper  conductors 
and  insulation  in  the  stator,  and  the  rotor.  The  total  "electromagnetic 
weight"  is  the  sum  of  the  weights  of  these  components.  Total  alternator 
weight  is  estimated  by  multiplying  the  electromagnetic  weight  by  1.20  to 
account  for  endbells,  bearings,  shafts,  and  structural  material. 

Magnetic  flux  lines 

Stator  conductors 


Insulation 


Figure  4.  Conventional  (Iron-Core)  Alternator  Radial  Cross  Section 
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Three  fundamental  relationships  are  used  in  sizing  alternators.  The 
rotor  diameter  is  determined  by  the  stress-limited  tip  speed,  Vt>  as 
fol lows , 

720  V. 

0  =  OPT  •  (inches>  (8) 

Machine  length  can  then  be  found  from  the  generated  phase  voltage  expres¬ 
sion  (Appendix  A), 

V  ,  =  7.4018  x  IQ'10  x  RPM  xK  xNxDxLxB,  (volts  rms )  (9) 

ph  w 

where  Kw  is  the  winding  factor,  N  is  the  number  of  series-connected  turns 
per  phase,  N  =  (p/PPATH)nsCs/2 ,  with  p  =  number  of  poles,  PPATH  =  number 
of  parallel  paths  per  phase,  ng  =  number  of  slots  per  pole  per  phase, 

Cs  =  number  of  conductors  per  slot,  L  is  the  stack  length  in  inches,  and 
B  is  the  peak  value  of  the  fundamental  sine  wave  component  of  flux  density 
in  the  air  gap  in  lines  per  square  inch.  The  phase  voltage  calculated 
from  Equation  9  is  the  open  circuit  phase  voltage.  Machine  reactance 
will  reduce  the  terminal  phase  voltage  under  load;  however,  since  machine 
reactances  are  not  calculated  in  the  algorithms,  this  voltage  regulation 
effect  has  not  been  included.  The  effect  of  machine  reactance  for 
sinusoidal,  steady-state  operation  is  to  reduce  the  phase  voltage  at  the 
terminals  to  a  value  given  by 

, . .  1/2 

\  -b  +  /  b2  -  U'(  p*/(9PF2)  | 

Vph  (LOAD)  "j  2  (  ’  (volts  rms)  (10) 


where  b  -  [ ( 2X&  P  /|_pp2  )/(3PF)]  -  V^,  Xg  is  the  synchronous 

reactance  in  ohms,  PF  is  the  power  factor  of  the  load,  P  is  the  output 
power  in  watts,  and  Vp^  is  given  in  Equation  9.  For  any  power  output 
greater  than  zero,  the  terminal  phase  voltage  under  load,  ^(load)* 
will  be  less  than  the  open  circuit  voltage  V  Thus,  machines  designed 
to  use  the  algorithms  presented  in  this  report  will  supply  the  rated 
power;  however,  the  rated  voltage  is  somewhat  high  and  the  rated  current 
is  somewhat  low,  due  to  the  fact  that  voltage  regulation  has  been  neglected. 
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For  a  given  power  output,  power  factor,  number  of  phases,  and  phase 
voltage  rating,  the  current  per  phase  in  the  stator  is  given  by 

Jph  =  P/Pf7^  Vph)  (amps  rms)  (11) 

A  generator  design  is  obtained  by  iterating  three  different, 
independent  design  variables.  These  are:  ng  (also  equal  to  the  number 
of  single  coils  in  a  pole  phase  group  as  indicated  in  Figure  5);  C&  (which 
is  always  an  even  number  for  a  double  layer  winding);  and  PPATH. 

The  number  of  slots  per  pole  per  phase,  ns,  is  varied  from  2  to  6 
with  ns  greater  than  or  equal  to  3  for  permanent  magnet  machines  to 
reduce  pulsation  losses.  Cs  is  varied  from  2  to  20  in  increments  of 
two.  PPATH  is  varied  from  1  to  n$p,  where  n$p  is  the  number  of  coils  per 
phase;  however,  only  those  values  which  divide  evenly  into  n$p  are 
considered.  For  example,  if  ngp  =  8,  only  the  values  1,  2,  4,  and  8  are 
used  for  PPATH.  A  double  layer  stator  winding  design  is  used  for  all 
iron-core  alternators,  as  illustrated  in  Figure  5.  Cg/2  conductors  lie 
in  an  upper  bar  and  C$/2  in  a  lower  bar  within  a  slot. 


[]  =  Phase  A  coil  side  Note:  ns  =  2  (two  slots,  or  coils,  per 
|  =  Phase  B  coil  side  pole  per  phase) 

Phase  C  coil  side 


electrical  degrees 


Figure  5.  Stator  Coil  Arrangement  for  a  Double  Layer  Winding 
With  Full  Pitch  Coils 
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Required  input  parameters  are  listed  in  Table  1.  Parameters  listed 
in  the  first  group  are  independent  input  variables  which  are  specified  to 
obtain  a  machine  of  the  desired  rating.  The  second  group  of  input  para¬ 
meters  are  set  by  the  present  technology  of  high  speed,  high-power- 
density  rotating  machinery.  For  example,  the  stator  bar  current  density 
is  limited  by  the  cooling  capability  of  the  stator.  Given  these  inputs, 
the  three  independent  design  variables  are  varied  to  find  a  minimum 
weight  machine.  Several  intermediate  parameters  of  interest  are: 

(a)  total  number  of  active  conductors  in  the  stator 

Z  =  ns  n^  p  C$  =  2  x  PPATH  x  n^  x  N  (12) 

(b)  total  number  of  slots 

SLOTS  =  n  nA  p  (13) 

s  <p 

(c)  current  per  conductor 

Ic  =  Iph/ppATH  (amps  rms)  (14) 

(d)  specific  electric  loading  (ampere  conductors  per  inch  of  stator 
periphery) 


0  =  _ £  =  .i_  _s.  _  x  _P _  x  (15) 

\  IID  PPATH  x  PFxV  .  x  720  V+  K  l 

ph  t 

(e)  slot  width  (Figure  3) 

s  =  HD  (a-1 )/(ansn^p) ,  (inches)  (16) 

where  at  is  the  ratio  of  the  magnetic  flux  density  in  the  tooth  to  that 
in  the  air  gap.  From  Figure  6(a),  it  can  be  seen  that  the  tooth  width 
is  given  by  s/(a-l)  when  the  assumption  is  made  that  the  slot  and  tooth 
have  straight  sides. 
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TABLE  1 

REQUIRED  INPUT  PARAMETERS  FOR  IRON-CORE  ALTERNATOR  DESIGNS 


PAREMETER 

SYMBOL 

UNITS 

NOMINAL  VALUE 

Rotor  Speed 

RPM 

rev. /min 

18,000 

Power  Rating 

P 

watts 

* 

Voltage  Rating  (phase) 

Vph 

volts  (rms) 

588 

Number  of  Phases 

"♦ 

- 

3 

Number  of  Poles 

P 

- 

14 

Power  Factor 

PF 

- 

* 

Rotor  Tip  Speed 

Vt 

ft/sec 

600 

Stator  Bar  Current  Density 

Jcu 

? 

amp/in 

20,000 

Peak  Gap  Flux  Density 

B 

2 

1 i nes/in 

50,000 

Insulation  Rating 

VPMIL 

vo Its /mil 

80 

Flux  Density  in  the  Yoke 

^yoke 

2 

lines/in 

120,000 

Stator  bar  Packing  Factor 

F 

- 

0.80 

Tooth  Flux  Density/B 

*Depends  upon  the  load. 

(f)  slot  depth  (see 

a 

Figure  6(a)) 

2.5 

C  I 

(s-26)  F  x  J 

—  +  0.01  (C 
cu 

s  -  2)  +  56, 

(inches) 

where  Jcu  is  the  current  density  in  the  copper  of  the  stator,  6  is 
the  ground  wall  insulation,  6  =  (Vp^/VPMIL)  x  1000,  and  F  is  the 
stator  conductor  packing  factor.  If  the  cross-sectional  area  of  a 
conductor  were  Ab>  then  the  copper  area  is  given  by  FAb  and  the  space 
reserved  for  cooling  passages  is  given  hy  ( 1  - F ) Ab . 


12 


AFAPL- TR- 79-2073 


h— s — 


t  =  s/(o-l) 


Figure  6(a).  Slot  Detail  with  Iron  Teeth 


Figure  6(b).  Stator  Coil  Geometry 


Figure  6.  Stator  Winding  Details  for  Iron-Core  Alternators 
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(g)  yoke  depth  (thickness  of  back-iron) 

dyoke  =  (D/p)  (B/Byoke^’  (inches)  08) 

where  B  is  the  allowable  flux  density  in  the  back-iron. 

(h)  end  turn  length  per  slot  (see  Figure  4) 

n  ?7r  1/2 

L'  -  21  -  -  x  [1  -  cos  ~  ]  ,  (inches)  (19) 

/?  sinx  P 

where  the  end  turn  extension  angle,  is  defined  in  Figure  6(b). 

(i)  stack  length  (straight  length  of  the  stator  coils) 

L  =  Vph/[7 . 401 8  x  10‘10  x  RPM  x  ^  x  (n^/p)  *  (p/PPATH)  x  D  x  B  x 

c  -0.06546 

^_Sj  (inches)  (20) 

where  the  factor  (CS/2)~^‘0B^B &  accounts  for  slot  leakage  as  the  number 
of  conductors  per  slot  (or  alternatively,  the  slot  depth)  is  increased. 
This  decreases  the  effective  flux,  B,  and  increases  the  required  machine 
length . 

(j)  frequency 

f  =  (p/2)({j{p)  •  (Hz)  (21) 

Once  the  input  parameters  of  Table  1  have  been  defined,  the  three 
independent  design  variables  n$,  Cs>  and  PPATH  are  varied  to  yield 
alternator  designs  which  consist  of  components  whose  weights  (in  pounds) 
are  calculated  by  the  following  formulas: 

(1)  yoke  (back-iron  weight) 

WV  ■  PY  t"d,oke  <°  +  2d  +  Ve)]  L'  (2Z> 
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(2)  teeth 

WT  -  PT  [  L,  (23) 

(3)  stator  insulation 

WINS  =  PINS  x  SL0TS  X  (sd  “  CsIc/Jcu)(L  +  2L,)’  (24) 

(A)  stator  copper 

“c  ■  »cu  *  SL0TS  *  <CsVJcu><L  *  L'>’  (25> 

(5)  rotor 

WR  =  pR  (ttD2/ 4 )  (1  -  S2)  L,  (26) 


where  the  mass  densities,  Py,py,  etc,  are  in  lbs/cu.  in.  Appendix  B, 
and  simply  multiply  the  calculated  component  volumes.  Overall  machine 
weight  is  given  By  (Wy+Wy+Wj^<.+Wcu+WR)  x  ^  and  the  "sPec’^c  weight" 
in  pounds  per  kilowatt  can  be  found  by  dividing  the  overall  machine 
weight  by  the  power  rating  in  kilowatts. 

An  algorithm  for  the  volume  of  the  generator  can  be  derived  once 
the  machine  dimensions  have  been  determined.  The  envelope  volume  in 
cubic  inches  can  be  found  by  multiplying  the  overall  machine  length 
by  the  frame  cross-sectional  area,  or 

VOL  =  ir  ( L  +  2e)  (D  +  2d  +  2dyQke  +  0.5)2/4,  (27) 

where  e  is  the  distance  which  the  end  turns  extend  beyond  the  straight 
portion  of  the  coil  (Figure  6(b)).  A  value  of  0.5  inch  is  assumed  to 
account  for  twice  the  air  gap  length  plus  twice  the  bore  seal  thickness. 

A  bore  seal  is  required  on  the  stator  inside  diameter  of  machines  which 
use  fluid-cooled  stators  in  order  to  exclude  coolant  from  the  air  gap. 

Appendix  B  includes  a  listing  and  sample  output  of  a  computer 
program  written  to  calculate  the  weight  and  volume  of  a  PMG.  Parametric 
results  obtained  from  the  program  are  described  in  Paragraph  3.4. 
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In  some  cases,  particularly  for  short  run  times,  it  may  be  desirable 
to  design  a  generator  without  active  cooling.  Generators  which  operate 
in  this  adiabatic  mode,  where  waste  heat  is  absorbed  by  the  heat  capacities 
of  the  materials,  are  called  "thermal  lag"  machines. 


Thermal  lag  designs  differ  from  actively  cooled  designs  due  to 
a  reduction  in  the  allowable  stator  conductor  current  density  and  an 
increase  in  F,  the  stator  bar  packing  factor.  For  an  operating  time 
of  A  t  seconds  and  an  average  temperature  rise  of  AT  (°F),  the  current 
density  can  be  obtained  from  the  following  equation: 


cu 


/ 


C,  AT  p 
h  cu 

At  p 

on 


(amp/i n* 


(28) 


where  C,  =  315  joules/lb/°F  is  the  heat  capacity  of  copper,  p  is 
n  ^  c  u 

0.32  lbs/in'. ,  and  p  is  the  electrical  resistivity  of  copper,  which  can 
be  calculated  from 


p  =  6.77  x  10'7 


»[ 


234.5  +  T 
254 . 5* 


(ohm-inches ) 


(29) 


where  T  is  in  degrees  centigrade.  The  average  temperature  rise  is 
assumed  to  be  about  1/4  of  the  allowable  hot  spot  (or  local)  temperature 
excursion,  which  may  be  almost  400°F;  therefore  AT  is  assumed  to  be 
100°F.  The  allowable  operating  stator  conductor  current  density  is 
only  about  8,500  amp/in  for  a  120  second  run  time.  Another  design 
consideration  for  thermal  lag  generators  involves  F.  For  actively  cooled 
designs,  about  20%  of  the  conductor  bar  cross  section  is  required  for 
cooling  passages,  while  0%  is  required  for  thermal  lag  designs.  Therefore, 
F  0.8  for  an  actively  cooled  generator,  while  F  =  1 .0  for  a  thermal 
machine . 
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3.3  GENERATOR  EFFICIENCY  CALCULATION 

An  estimate  of  generator  efficiency  is  obtained  by  calculating  the 
losses  due  to  ohmic  heating  in  the  stator  conductors  and  eddy  current 
and  hysteresis  losses  in  the  iron  stator  yoke  and  teeth.  This  approxi¬ 
mation  should  be  accurate  to  within  10X,  since  the  rotor  losses, 
including  induced  surface  current  losses,  bearing  and  windage  losses, 
are  much  less  than  the  stator  losses  due  to  the  large  operating  current 
densities  in  the  stator  bars,  high  operating  frequency  and  relatively 
high  magnetic  flux  densities  in  the  iron. 

Ohmic  heating  in  the  stator  windings  can  be  calculated  from 

Pohi»ic  ■  Jcu  »  «o  •  <"“«>  (30> 

2 

where  Jcij  is  the  copper  current  density  in  amp/in  ,  p  is  the  electrical 
resistivity  of  copper  defined  in  Equation  29  and  Vcu  =  Wcu/Pcu  is  the 
volume  of  stator  copper  in  cubic  inches. 

Iron  losses  consist  of  hysteresis  losses  and  eddy  current  losses 
given  by 

Ph  =  7.1  x  10'5  f  B1 -6  (watts/lb)  (31 ) 

and 

Pc  =  4.3  x  10'9  f2  B2  (watts/lb)  (32) 


where  f  is  the  operating  frequency  and  B  is  the  peak  magnetic  flux 
density.  Actual  loss  in  the  yoke  or  teeth  can  be  calculated  by  multiplying 
the  loss  per  pound  by  the  weight  of  the  yoke  or  teeth  in  pounds. 

Generator  efficiency  is  given  by 


nG 


out 
+  P 


loss 


(33) 


where  P  t  is  the  generator  output  power  in  watts  and  P-joss  is  the 
sum  of  the  ohmic  losses  and  the  iron  losses. 
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3.4  COMPUTER  PROGRAM  PARAMETRIC  RESULTS 

To  determine  the  accuracy  of  the  computer  program  which  utilizes 
the  equations  developed  in  Paragraph  3.2,  the  results  obtained  from  the 
computer  program  are  compared  in  Table  2  with  an  actual  high  power  PMG 
design  performed  by  Ai Research  Manufacturing  Company  of  California.  The 
calculated  specific  weight  of  0.104  Ib/KW  agrees  with  the  design  value, 
and  the  calculated  volume  is  very  close  to  the  design  value.  Of  course, 
all  of  the  independent  triable  values  and  design  parameters  were  chosen 
to  be  the  same  for  both  the  computer  program  and  the  Ai Research  design. 
This  design  can  be  considered  as  a  "base-line  case"  with  the  parametric 
results  described  below  obtained  by  varying  the  independent  variables  P 
(power),  Vph  (phase  voltage),  or  RPM  (rotational  speed)  one  at  a  time 
while  keeping  the  other  variables  equal  to  their  base  value.  In  addition, 
parameters  which  depend  upon  the  state-of-the-art,  including  Jcu  (stator 
conductor  current  density)  and  (rotor  tip  speed),  are  varied,  in  turn, 
from  their  base  values  while  keeping  the  other  variables  at  their  base 
value. 


Before  describing  the  parametric  results  obtained  from  the  computer 
program  it  is  necessary  to  show  that  a  fair  comparison  of  generator 
weights  and  volumes  for  different  values  of  the  independent  variables 
can  be  made  only  if  the  different  designs  have  the  same  value  of  specific 
electric  loading,  in  ampere-conductors  per  inch  of  stator  periphery, 
given  by 


n  .  C  *  p 
s  s  r 

\  =  720- PPATH -PF 


'P  (RPM)  ‘ 

Vnh -V. 

L  Ph  t  J 


(34) 


Since  the  number  of  poles,  p,  and  the  load  power  factor,  PF,  are  not 
varied  in  this  study,  the  design  parameters  n$ ,  Cs ,  and  PPATH  must  vary 
to  yield  designs  with  comparable  specific  loadings  for  different  values 
of  the  independent  variables  in  brackets.  For  example,  if  P  were  in¬ 
creased  from  10  to  20  MW,  the  loading  would  double;  however,  if  n$  C  / 
PPATH  were  correspondingly  halved,  the  loading  would  remain  the  same 
and  a  fair  comparison  of  generator  specific  weight  and  volume  for  designs 
at  the  two  power  levels  could  then  be  made. 
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TABLE  2 

BASE-LINE  COMPARISON,  5  MW  AIRESEARCH  PERMANENT  MAGNET  GENERATOR 


Design  Parameters 

Ai Research  Design 

PMG  Program 

Voltage  (Volts,  rms,  line- line) 

1018.5 

1018.5 

Power  (MW) 

5 

5 

RPM 

18000.0 

18000.0 

No.  of  phases 

3 

3 

No.  of  poles 

14 

14 

Rotor  tip  speed  (ft/sec) 

625.50 

625.50 

Power  factor 

0.82 

0.82 

Armature  current  density 

36,270 

36,270 

(amp/in^) 

2 

Peak  gap  density  (lines/in  ) 

48,900 

47,400 

Insulation  rating  (wolts/mil) 

100 

100 

Yoke  flux  density  (lines/in  ) 

100,000 

100,000 

Winding  pitch 

5/6 

5/6 

Ratio  of  tooth  flux  density  to 

2.2 

2.2 

air  gap  flux  density 

No.  of  conductors  per  slot 

2 

2 

No.  slots/pole/phase 

4 

4 

No.  of  slots  in  the  stator 

168 

168 

No.  of  parallel  paths 

14 

14 

Specific  electric  loading 

3177 

3120 

(ampere-conductors/inch ) 
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TABLE  2  (CONTINUED) 


Weiqhts  (lbs) 

Ai Research  Design 

PMG  Program 

Total  machine 

520 

522 

Total  machine  (Ib/kw) 

0.104 

0.104 

Stator  copper 

26.2 

24.8 

Back  iron  (yoke) 

64.6 

65.5 

Rotor 

309.4 

308.7 

Dimensions  (inches) 

AiResearch  Desiqn 

PMG  Proqram 

Stator  coil  end  extensior 

i  1  .01 

1  .06 

Stack  length 

31  .22 

31  .19 

Stack  OD 

9.515 

9.08 

Rotor  OD 

7.96 

7.96 

Volume  (cubic  feet) 

2.18 

2.03 

Aspect  ratio  (L/D) 

3.92 

3.92 

Losses  (Kilowatts) 

AiResearch  Desiqn 

PMG  Proqram 

Iron  losses 

Back-iron  (yoke) 

27.8 

27.9 

Teeth 

18.2 

13.7 

Ohmic  losses 

135.8 

126.8 

Efficiency  {%) ,  based 

on  stator  losses 

only 

96.5 

96.7 

20 
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Specific  electric  loading  is  an  indication  of  how  hard  a  machine 
is  being  worked  (or  loaded).  Higher  specific  loadings  imply  more  power 
output  per  pound  of  machine  weight  as  shown  in  Figure  7.  The  loading 
cannot  be  increased  indefinitely  because  of  cooling  considerations.  For 
advanced  high-power  permanent  magnet  generators,  the  specific  electric 
loading  is  less  than  3500.  Note  also  in  Figure  7  that  for  each  value 
of  loading,  i.e.,  1560,  2340,  or  3120,  the  specific  weight  depends  upon 
the  particular  combination  of  the  two  design  parameters  n$ ,  and  C  .  The 
dependence  is  not  strong,  therefore  n  and  C  can  be  varied  (to  obtain 

5  S 

the  same  value  of  loading)  as  the  independent  variables  are  changed  and 
a  fair  comparison  of  generator  weights  and  volumes  can  still  be  made. 

In  Figure  8(a),  the  specific  weight  and  volume  for  high-power  per¬ 
manent  magnet  generator  designs  are  shown  as  a  function  of  rated  output 
power.  Note  that  the  specific  weight  decreases  even  out  to  20  MW; 
however,  the  decrease  is  much  less  for  power  levels  above  5  MW.  It  is 
doubtful  that  generators  above  10  MW  can  be  constructed  as  single  units 
because  of  manufacturing  and  assembly  problems  and  extremely  large  aspect 
ratios.  An  approximate  value  for  the  aspect  ratio  (stack  length 
divided  by  rotor  OD)  can  be  obtained  from  the  following  equation. 

L/D  :  70  -  p-  (-Rpm2) -  (35) 

\  \  B  Ql  PF 

Note  from  Figure  8(b)  that  the  aspect  ratio  is  approximately  4:1  at 
5  MW  and  increases  linearly  with  power,  P,  with  the  other  parameters 
fixed  at  the  base-line  values.  One  way  to  reduce  the  aspect  ratio  is 
to  increase  the  rotor  tip  speed;  however,  AiResearch  has  chosen  625.5  feet 
per  second  as  a  fairly  conservative  number  which  has  been  experimentally 
demonstrated.  The  approach  taken  by  AiResearch  is  to  design  5  MW 
machines  as  modules  which  are  connected  together  in  order  to  achieve 
higher  power  ratings.  The  same  approach  will  be  taken  in  this  report 
with  a  20%  weight  penalty  and  a  50%  volume  penalty  imposed  for  inter¬ 
connection  as  described  in  Section  V.  Note  also  in  Figure  8(a)  the 
very  small  volume  of  the  PMG.  A  5  MW  machine  occupies  only  2  cubic 
feet  of  space. 
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Figure  8(a).  Specific  Weight  and  Volume  Versus  Output  Power 


P  (Output  Power  in  Megawatts) 

Figure  8(b).  Aspect  Ratio  (L/D)  Versus  Output  Power 
Figure  8.  PMG  Output  Power  Variation 
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In  Figure  9,  the  specific  weight  and  volume  are  shown  as  a  function 
of  the  generator  1 i ne-to-neutral  voltage  rating.  Note  the  weak  dependence 
below  about  1  KV  and  the  significant  upswing  above  this  rating.  Designs 
above  about  1.2  KV  were  impractical  because  there  was  no  room  for  copper 
left  in  the  slot  due  to  the  thickness  of  insulation  required.  In 
Figure  10,  the  specific  weight  and  volume  are  shown  to  decrease  with 
larger  rotational  speeds.  The  effect  on  specific  weight  and  efficiency 
is  indicated  in  Figure  11  for  two  different  rotor  tip  speeds.  Increased 
tip  speed  results  in  lighter  and  more  efficient  generators.  If  the 
operating  current  density  in  the  stator  conductors  is  varied,  the 
specific  weight,  volume,  and  efficiency  are  affected  as  shown  in  Figure  12. 
In  general,  higher  current  densities  result  in  smaller  and  lighter 
generators;  however,  the  efficiency  is  reduced.  The  efficiency  decreases 
because  the  conduction  (I  R)  losses  increase. 


Figure  12.  PMG  Stator  Current  Density  Variation 
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SECTION  IV 

SUPERCONDUCTING  HIGH-POWER  GENERATORS 


4.1  BACKGROUND 

The  equations  developed  are  restricted  to  that  class  of  super¬ 
conducting,  synchronous  AC  generators  which  appear  to  be  lightest  at 
power  levels  around  20  MW,  i.e.,  image-shielded  machines  which  use  a 
conductive  environmental  shield  to  attenuate  the  rotating  magnetic  field 
in  the  aircraft  environment  and  an  actively  cooled  stator  operating  at 
high  current  densities.  This  kind  of  superconducting  generator  is  being 
developed  under  contract  to  General  Electric  Company  as  a  second 
generation  machine  which  utilizes  advanced  concepts  and  technology  to 
operate  within  requirements  such  as  rapid  start-up  (one  second  excitation 
and  spin-up)  and  rectified  output  (References  7,  8). 

4.2  WEIGHT  AND  VOLUME  EQUATIONS 

A  radial  cross  section  of  the  superconducting  generator  is  shown  in 
Figure  13.  The  distance  from  the  superconducting  coil's  outside  radius 
to  the  stator  conductors  can  be  considered  as  the  "air  gap,"  since  no 
iron  is  present.  The  pertinent  tip  speed  is  the  peripheral  speed  at  the 
winding  outer  radius.  This  radius  is  given  by 

720  x  \J 

Rr  =  ~2tt~R'PM  (inches)  (36) 

where  Vt  is  the  tip  speed  (ft/sec)  and  RPM  is  the  operating  rotational 
speed . 

A  longitudinal  cross  section  of  the  machine  is  shown  in  Figure  14, 
where  all  of  the  relevant  dimensions  are  indicated.  The  stator  inside 
diameter,  D,  which  is  used  below,  can  be  calculated  from 

D  =  2  [Rr  +  At  +  As  +  g]  (inches)  (37) 
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Figure  13.  Superconducting  Generator  Radial  Cross  Section 


29 


AFAPL-TR-79-2073 


In  this  case,  D  cannot  be  assumed  to  be  the  same  as  D  ,  i.e.,  2R  , 

r  r 

because  of  the  large  "air  gap."  It  should  be  noted  that  there  are 
no  high  permeability  materials  within  the  flux  paths  of  this  kind  of 
machine;  however,  the  superconducting  field  windings  can  easily  supply 
the  large  mmf  required  by  this  low  permeability  path  by  virtue  of  the 

extremely  large  current  densities  in  the  superconducting  field  winding 

2 

modules  (96,774  amp/in  ). 

The  electromagnetic  length  of  the  stator  can  be  obtained  by  using 
the  voltage  generation  equation  as  in  the  previous  section. 

Lm  =  P/[7 .50  x  10'6  x  PF  x  Kw  x  RPM  x  D2  x  Bkg  x  Qj  (inches)  (38) 

The  term  B^g  is  the  average  flux  density  (crest  value)  in  Kilogauss  in 
the  stator  winding  calculated  by  taking  an  average  of  the  two  crest 
values  at  the  stator  winding  ID  and  OD.  This  must  be  done  since  there 
are  no  iron  teeth  between  stator  bars  in  a  superconducting  generator, 
hence  the  permeability  of  the  magnetic  circuit  is  low  and  the  magnetic 
flux  density  decreases  rapidly  across  the  winding.  A  subroutine  is  in¬ 
cluded  as  a  part  of  the  weight  and  volume  algorithm  which  calculates  flux 
densities  based  on  the  Biot-Savart  law.  Only  four-pole,  rectangular 
cross-sectional  field  coils  can  be  handled  at  the  present  time.  Specific 
electric  loading  Q^,  is  the  same  term  used  in  Section  III;  however,  it  is 
more  convenient  in  the  present  context  to  express  as 

Ql  =  AaF'Jcu  (amp  conductors  per  inch)  (39) 

where  is  the  winding  thickness,  F1  is  the  armature  (or  stator)  packing 
factor  and  J  is  the  copper  current  density.  F 1 J  is  the  armature 
overall  annulus  current  density.  In  the  computer  program  in  Appendix  B, 

F1  is  represented  by  the  variable  FARM  and  is  a  calculated  quantity. 

The  required  length  of  the  active  conductor  straight  section  is  less  than 
Lm  because  of  emf  generated  in  the  end  turns  for  this  kind  of  machine. 

L,  the  actual  straight  section  length,  is  given  by  the  following  equation, 
where  Lm  is  assumed  to  be  25%  longer  than  the  actual  length. 

L  =  Lm/1.25  (inches) 


(40) 
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The  stator  coil  geometry  is  identical  to  that  in  Figure  6(b)  and  the 

same  geometrical  relationships  apply.  However,  a  single  layer  stator 

winding  is  assumed  as  shown  in  Figure  15(a).  A  detail  of  a  single  bar 

(which  occupies  a  single  slot  on  the  stator)  is  shown  in  Figure  15(b). 

The  ratio  of  cross-sectional  area  of  copper  to  the  bar  cross  sectional 

area  is  the  bar  packing  factor  F^ar-  A  hollow,  water-cooled  conductor 

is  indicated  and  the  cross-sectional  area  of  copper  is  just  C$ Ic/^cu , 

where  C  is  the  number  of  conductors  per  bar  and  I  /J  is  the  cross- 
s  r  c  cu 

sectional  area  of  a  single  conductor  carrying  a  current  of  I  amps  with 

2  c 

a  current  density  in  the  copper  of  JC(j  amps/in  . 

Racetrack-shaped  coil  modules  are  assumed  for  the  superconducting 
field  windings  as  indicated  in  Figure  16.  Module  dimensions  depend 
upon  Dr>  p,  and  L.  The  independent  dimensions  of  the  module  are  its 
height,  a,  and  bend  radius,  b.  As  shown  in  Figure  17,  once  a  and  b  are 
specified,  the  module  width  can  be  calculated  from 

t  =  [(R  -  a  cos  it/ p )  x  sin  n/p]  -  b  (inches)  (41) 

For  the  image-shielded  machine,  the  skin  depth  in  the  environmental 
shield  material  at  the  operating  frequency  must  be  calculated.  The 
frequency  is  given  by  Equation  21,  and  the  skin  depth  by 

5  =  0154  x>/^%  (inches)  (42) 

*  o 

where  u  =  2irf,  p  =  4u  x  10"^  and  o  is  the  conductivity  of  the  shield 
0  7 

material,  which  is  about  3.77  x  10  mhos/m  for  aluminum.  The 
environmental  shield  thickness,  Ags,  is  assumed  to  be  three  times  the 
skin  depth  for  adequate  attenuation  of  the  electromagnetic  field  within 
the  aircraft  environment. 

I 

i 


AFAPL-TR- 79-2073 


Figure  15(a).  Single  Layer  Stator  Winding  (Developed  View) 


Figure  15(b).  Stator  Slot  Detail 
Figure  15.  Stator  Winding  Details  for  an  Air-Core  Alternator 
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The  required  input  parameters  are  listed  in  Table  3.  Intermediate 
parameters  required  to  calculate  component  volumes  are  listed  below  with 
all  dimensions  in  inches,  unless  otherwise  indicated. 


> 


VDC/2-3 


(volts  rms) 


‘Ph  *  wf  %  V 

Rb  =  (Rr  "  a  cos  tt/P )/2 


rs  =  Rb  +  Ab 


R  =  +  At  +  A  +  g  +  A  +  s 


he  =  Rs  '  (s  +  rs) 


Lc  =  L  +  2  (b  +  t) 


L  +  2e 


SW(slot  width)  =  [ tt D  -  n^p  -  (ns'l  )nl+,P(5t)ar]/SL0TS 


Aa  =  CsIph/(SW-261)  +  °’01  (Cs_1)  +  26i 


Lde  =  V5 


L  ,  =  L  /5 

ade  c' 


z  -  mL/iph 


Acu  =  (amp/i"2) 


(43: 

(44; 


(45; 

(46: 

(47! 

(48) 

(49; 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


There  are  thirteen  components  for  which  volumes  are  calculated. 
Weights  of  the  individual  components  are  then  determined  by  multiplying 
the  component  volumes  by  the  corresponding  mass  densities  shown  in 
Table  8  in  Appendix  B. 
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TABLE  3 

INPUT  PARAMETERS  REQUIRED  FOR  CALCULATING  THE  WEIGHT 


AND  VOLUME  OF  A 

SUPERCONDUCTING  GENERATOR 

PARAMETER 

SYMBOL 

UNITS 

NOMINAL  VALUE 

Power  rating 

P 

Watts 

20  x  106 

Field  winding  tip  speed 

Vt 

ft/sec 

471 

RPM 

RPM 

revolutions/min 

6,000 

Stator  copper  current  density 

Jcu 

2 

amp/in 

33,000 

Field  winding  current  density 

Jf 

2 

amp/cm 

15,000 

DC  load  voltage 

VDC 

volts 

40,000 

Power  factor 

PF 

-- 

0.86 

Winding  distribution  factor 

-- 

1  .00 

Stator  bar  packing  factor 

Fbar 

— 

0.45 

Number  of  poles 

P 

— 

4 

Number  of  phases 

0, 

<f> 

-- 

3 

Coil  bend  radius 

b 

inches 

1  .50 

Coil  module  height 

a 

inches 

2.85 

Stator  coil  end  turn  angle 

X 

degrees 

45 

Gap  between  rotor  and  stator 

g 

i nches 

0.65 

Distance  from  stator  OD  to 
environmental  shield  ID 

s 

inches 

6.0 

Damper  shield  thickness 

As 

inches 

0.75 

Phase  break  insulation 

6ph 

inches 

0.50 

Stator  bar  insulation 

6bar 

inches 

0.02 

Stator  conductor  insulation 

6i 

i nches 

0.005 

Bore  tube  thickness 

Ab 

inches 

0.25 

Torque  tube  thickness 

At 

i nches 

0.60 

End  bell  thickness 

Ae 

inches 

0.50 
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Component  volumes,  in  cubic  inches,  are  given  by  the  following  equations: 

(1)  Environmental  shield  —  cylindrical  portion 

V1  =  "Lc  x  {2Aes  Rs  +  Aes}  (57) 

(2)  Environmental  shield  -  two  conical  sections 

V2  =  2ttS  Aes  x  (2Rs  +  Aes  '  s)  (58) 


(3)  End  Bells 

V3  =  2tt  A0  x  [(Rs  -  s)2  -  r2] 

(4)  Superconducting  Winding  Modules 
V4  =  p  x  [2Lat  +  ua  (2bt  +  t2)] 

(5)  Drive  end-stub  shaft  plus  bearings 


(59) 


(60) 


(61) 


(6)  Anti-drive  end  —  helium  transfer  system  plus  bearings 

V6  =  7rrsLade  (62) 

(7)  Torque  tube  plus  torque  tube  extension 

V7  =  tt(Lc  +  s)  (2Rr  At  +  l\)  (63) 

(8)  Bore  tube 

V8  =  wLc  (2Rb  Ab  +  Ab>  (64) 


(9)  Field  winding  support  structure 
V9  =  *Lc  (Rr  -  rs}  '  V4 


(65) 


(10)  Electromagnetic  shield  on  the  rotor  (damper  shield) 

let  R  =  R  +  A*  and  R  =  R  +  A  , 
es  r  t  eso  es  s 


V’"  ■  ’Lc  '“eso  -  Res>  +  *4s<4o  '  's* 


'10 


(66) 


38 


AFAPL- TR- 79-2073 


(11)  Copper  in  the  stator 

Vn  =  ZACU  (L  +  W/sinx)  (67) 

(12)  Insulation  in  the  stator 

let  Ri  =  Reso  +  g*  and  Ro  =  Ri  +  V 

V12  =  ttL'  (R2  -  R2)  .  vn  (68) 

(13)  Bore  seal 

V13  =  1tL'  ^2dBS  x  Ri  +  dBS^’ 

where  dg,.  is  the  bore  seal  thickness  and  is  calculated  in  the  program. 

A  computer  program  was  written  to  calculate  the  volumes  and  the 
corresponding  weights.  The  total  weight  of  the  machine  is  just  the  sum 
of  the  individual  component  weights.  No  correction  factor  is  required 
to  account  for  end  bells,  bearings,  structure,  etc.,  as  in  the  PMG 
design,  since  everything  is  included  in  the  component  weights.  An  enve¬ 
lope  volume  which  includes  the  cylindrical  portion  of  the  environmental 
image  shield  and  the  two  end  pieces  which  are  frustums  of  cones  is 
calculated  from 

VOLUME  =  7tL  1  (R  +  A  )2  +  27rS[2Rf  +  S2  -  2SR  +  R  (R  -S)].  (cu/in)  (70) 

5  GS  5  5  5  S 

Appendix  B  includes  a  listing  of  the  computer  program  and  a  sample  output. 
Parametric  results  obtained  from  the  program  are  described  in  Paragraph  4.4 

4.3  GENERATOR  EFFICIENCY  CALCULATION 

An  estimate  of  the  efficiency  for  the  superconducting  generator 
can  be  obtained  by  calculating  the  ohmic  (I  R)  losses  in  the  normally 
conducting  copper  stator  bars,  the  eddy  current  losses  in  the  stator 
bars,  and  the  losses  due  to  surface  currents  induced  on  the  environmental 
image  shield.  Even  though  windage,  bearing,  and  rotor  (electromagnetic 
shield)  losses  are  not  considered,  they  are  small  compared  to  the  above 
losses,  and  a  reasonably  accurate  estimate  of  efficiency  should  be  obtained 
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Ohmic  losses  are  determined  by 


P  ..  =  iT  p  V-,-,  (watts) 

ohmic  cu  c  11 


where  Jcu  is  the  current  density  in  the  stator  in  amperes/m  ,  pc  is  the 
resistivity  of  the  conductors,  assumed  to  be  copper,  given  by 

Pc  =  1.72  x  10  8  [^254  5~^]  (ohm-meters)  [72 


where  T  is  the  average  stator  bar  conductor  temperature  in  degrees 
centigrade  and  is  the  volume  of  copper  in  the  stator  given  by 
Equation  67,  expressed  in  m3. 

Since  the  full  magnetic  field  is  allowed  to  penetrate  the  copper 
stator  bars  rather  than  being  shunted  through  iron  teeth  around  the  bars 
as  in  conventional,  iron-core  machines,  the  eddy  current  losses  are 
significant.  For  square,  hollow  conductors  as  shown  in  Figure  15b,  the 
eddy  current  loss  in  watts  is  given  by  (Reference  12) 


2  [V-h4 

to  o  l 


eddy  24p  .2  _  ,2 

o  i 


x  [B3  +  B3]  x  Z  x  A  x  L 

L  0m  rmJ  cu  c 


where  ta  is  the  electrical  radian  frequency,  hQ  is  the  conductor  height 

in  meters,  t^  is  the  width  of  the  hollow  center  of  the  conductor  in 

meters,  8&m  is  the  average  of  the  peak  value  of  the  azimuthal  components 

of  magnetic  flux  density  calculated  at  the  inner  and  outer  radii  of  the 

stator,  and  Brm  is  the  average  of  the  radial  component.  BQm  and  Brm 

are  in  Tesla.,  and  A  and  L  must  be  in  meters. 

V  cu  c 

Losses  due  to  induced  surface  currents  on  the  environmental  shield 
are  calculated  from 


shield 


x  ,  (watts) 
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where  R$  is  the  radius  of  the  environmental  shield  (in  meters),  L‘  =  L+2e 
(in  meters)  and  Btan  is  the  crest  value  of  the  azimuthal  component  of  the 
magnetic  flux  density  at  radius  Rg  in  Tesla.  The  skin  depth,  <5,  must  be 
expressed  in  meters.  Shield  losses  predicted  from  this  equation  are 
almost  50%  lower  than  the  losses  predicted  by  the  GE  design  as  shown 
in  Table  4. 


Generator  efficiency  is  calculated  from 


n 


G 


P  . 
out 


P  P 
out  +  loss 


(75) 


where  PQUt  is  the  electrical  power  output  and  P-|oss  is  equal  to  the  sum 
of  the  stator  losses  and  the  losses  in  the  environmental  shield. 

4.4  COMPUTER  PROGRAM  PARAMETRIC  RESULTS 

To  determine  the  accuracy  of  the  computer  program  which  implements 
the  weight  and  volume  equations  given  in  Paragraph  4.2,  predicted  machine 
parameters  from  the  program  are  compared  in  Table  4  with  an  actual 
superconducting  generator  design  produced  by  General  Electric  Company 
(Reference  8).  As  shown  in  Table  4,  the  predicted  weight  and  volume 
agree  well  with  the  GE  design  numbers  for  identical  input  parameters. 

This  design  point  is  considered  the  "base-line  design"  and  parametric 
results  are  obtained  by  varying  certain  input  parameters  from  their  base 
values,  indicated  in  Table  4. 

Parametric  results  are  obtained  for  the  independent  variables  P 
(output  power  rating),  DC  rectified  output  voltage,  and  RPM.  In  addition, 
the  technological  parameters  Jcu  (stator  conductor  current  density)  and 
(field  winding  current  density)  are  varied  from  the  base-line  design. 

As  described  in  the  beginning  of  Paragraph  3.4,  all  designs  must  be 
compared  on  an  equal  specific  electric  loading  basis.  This  has  been  done 
in  the  following  curves  which  illustrate  how  the  weight,  volume,  and 
efficiency  of  this  class  of  superconducting  generators  change  with  rating 
and  other  machine  input  parameters. 
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TABLE  4  (CONTINUED) 


Weiqhts  (lbs) 

GE  Design 

SCGEN  Program 

Total  machine 

1740 

1706 

Total  machine  (lb/KW) 

0.087 

0.085 

Rotor 

870 

846 

Armature  copper 

286 

274 

Environmental  shield 

280 

285 

Stator  bore  seal 

32 

39 

EM  shield 

112 

85 

Field  winding  modules 

246 

258 

Dimensions  (in.) 

GE  Desiqn 

SCGEN  Program 

Field  winding  outer  radius 

9.0 

9.0 

Armature  thickness 

1 .0 

1 .42 

Straight  length  of  coils 

6.0 

5.0 

Thickness  of  environmental  shield 

0.50 

0.69 

Radius  of  environmental  shield 

19.7 

19.1 

Volume  (cubic  feet) 

35 

35 

Losses  (KW) 

GE  Design 

SCGEN  Program 

Ohmic  losses 

770 

756 

Armature  eddy  current 

390 

347 

Environmental  shield 

380 

193 

Efficiency  {%) ,  based  on 
stator  losses  only 

93 

93.9 

Magnetic  fields  (kilogauss) 

GE  Design 

SCGEN  Program 

Radial  flux  density  at: 

armature  inner  radius 

14.0 

17.3 

armature  outer  radius 

9.9 

11 .6 
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In  Figure  18,  the  specific  weight,  volume,  and  efficiency  are  shown 
as  a  function  of  rated  output  power.  The  efficiency  of  this  class  of 
generators  is  comparable  to  the  PMG  only  at  power  levels  above  50  MW. 

Also,  the  volume  is  approximately  three  times  as  great  as  the  volume  of  a 
PMG  of  similar  power  rating. 

The  dependence  of  specific  weight,  volume,  and  efficiency  upon  the 
generator  output  voltage  is  indicated  in  Figure  19.  Phase  voltage  (rms) 
is  approximately  one  half  the  DC  rectified  output  voltage.  Note  that 
even  for  an  output  voltage  of  100  K V  DC  the  specific  weight  increases  by 
only  50%  above  the  base-line  design  at  40KV,  and  the  volume  by  only  27%. 
This  is  in  direct  contrast  to  the  PMG  where  voltages  of  even  a  few  KV 
were  unattainable.  This  high  voltage  capability  of  the  superconducting 
generator  is  due  to  the  "air-core"  nature  of  the  stator  winding  and 
consequently  the  extra  space  for  insulation  in  the  absence  of  iron  teeth. 

It  is  concluded,  therefore,  that  the  PMG  is  a  "low  voltage"  generator 
while  the  superconducting  generator  is  a  "high  voltage"  generator. 
Efficiency  is  only  a  weak  function  of  output  voltage  as  shown  in  Figure  19. 

The  effect  of  rotational  speed  variation  is  shown  in  Figure  20.  Note 
that  there  is  a  minimum  in  the  specific  weight  and  volume  curves  at 
about  8000  RPM  and  a  peak  in  the  efficiency  curve  at  about  the  same  speed. 

The  current  density  in  the  field  winding  coils  is  determined  by  the 
characteristics  and  packing  factor  of  the  superconducting  wire  or  cable 
used  in  winding  the  coils.  This  parameter  is,  therefore,  set  by  current 
technology  and/or  manufacturing  limits.  For  the  base-line  design,  the 
current  density  is  15,000  amps/cm  (i.e.,  96,774  amp/in  ).  In  Figure  21, 
the  generator  specific  weight,  volume  and  efficiency  are  shown  as  a 
function  of  field  winding  current  density.  Specific  weight  and  volume 
decrease  slowly  as  the  field  winding  current  density  is  increased.  Note 
that  the  efficiency  also  decreases.  This  offsets  the  advantage  of 
reduced  weight  and  volume  to  some  degree.  The  decrease  in  efficiency 
is  due  to  increased  eddy  current  losses  in  the  stator  winding  and 
environmental  shield. 


44 


Specific  Weight  (lb/KW) 


AFAPL-TR-79-2073 


Base-line  design  parameters : 
6000  RPM 


17,378  volts/phase 


_97 


_96 

_95 


_94 

Efficiency 

(*) 

_93 


1—92 


Ml 


90 


Figure  18.  Superconducting  Generator  Output  Power  Variation 
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In  Figure  22,  generator  specific  weight,  volume,  and  efficiency  are 
plotted  as  a  function  of  the  current  density  in  the  stator  conductors. 
Again,  all  three  parameters  decrease  for  increasing  current  density; 
however,  the  rate  of  decrease  of  both  weight  and  volume  is  greater  than 

that  observed  when  the  field  winding  current  density  was  varied. 

2 

Decreasing  efficiency  is  caused  by  the  larger  conduction  (I  R)  losses  in 
the  stator  conductors. 
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SECTION  V 

GENERATOR  INFLUENCE  ON  POWER  SYSTEM  WEIGHT 

A  true  comparison  of  different  generator  concepts,  such  as 
permanent  magnet  and  superconducting,  can  only  be  made  by  comparing 
complete  power  supply  systems.  Even  though  one  generator  may  weigh 
less  than  another  generator  for  the  same  power  output,  the  overall  system 
which  uses  the  lighter  generator  as  a  component  may  be  appreciably 
heavier  because  of  the  influence  that  the  generator  has  upon  the  system. 
The  mechanical  part  of  the  system  is  influenced  by  the  generator  (rotor) 
rotational  speed  and  the  generator  efficiency.  The  cooling  subsystem  is 
also  affected  by  generator  efficiency.  Components  in  the  electrical  part 
of  the  system  (power  conditioning  subsystem)  are  influenced  by  the 
generator  frequency  and  output  voltage. 

Therefore,  to  compare  the  PMG  and  the  superconducting  generator, 
system  weights  were  calculated  by  using  the  algorithms  described  in 
Appendix  C  (Reference  13)  and  the  two  generator  algorithms  in  Appendix  B. 
The  systems  are  shown  in  block  diagram  form  in  Figure  23  for  the  PMG  and 
Figure  24  for  the  superconducting  generator.  A  total  power  generation 
time  (or  run  time)  of  120  seconds  is  assumed.  The  auxiliary  mechanical 
power  is  assumed  to  be  zero  for  simplicity  and  a  main  transformer  (Tl) 
is  not  required  for  the  system  in  Figure  24  because  of  the  high  voltage 
capability  of  the  superconducting  generator.  The  assumed  component 
efficiencies  are  indicated  below  each  component.  Note  that  the  efficiency 
for  the  power  conditioning  component  includes  the  efficiencies  of  both 
a  rectifier  and  an  LC  filter.  Power  levels  shown  in  parentheses  are  for 
a  20  MW  (approximately)  generator  and  the  other  power  levels  correspond 
to  a  10  MW  (approximately)  generator.  The  load  consists  of  a  main  power 
of  6  MW  or  15  MW  at  30  KV  DC  or  40  KV  DC,  respectively,  and  a  secondary 
power  of  1  MW  or  2  MW  at  130  KV  DC  or  150  KV  DC,  respecti vely .  Numbers 
in  circles  near  the  gear  box  indicate  point  designs  which  use  a  gear  box 
or  direct  drive  between  the  turbine  and  the  alternator. 
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A  summary  of  the  seven  point  designs  is  shown  in  Table  5  along 
with  the  two  base-line  designs.  Generator  power  requirements  are 
determined  by  the  load  and  the  component  efficiencies  between  the 
generator  and  the  load.  The  superconducting  generator  output  power 
of  8.670  MW  is  slightly  less  than  the  PMG  output  power  requirement 
of  9.066  MW  (two  4.533  MW  modules),  since  the  inefficiency  of  a  high 
power  transformer  is  absent  in  the  case  of  the  superconducting  generator 
system.  The  10  MW  designs  include  a  directly  driven  and  a  gear  box 
driven  PMG  and  a  gear  box  driven  superconducting  generator.  Note  that 
two  PMG's  must  be  used  to  meet  the  power  requirement  since  it  has  been 
assumed  that  single  unit  lightweight  permanent  magnet  generators  larger 
than  about  5  MW  are  not  presently  practical;  however,  the  combination  of 
any  number  of  these  size  devices  (with  a  20%  weight  and  50%  volume 
penalty  for  connection)  is  feasible.  The  20  MW  point  designs  include 
a  direct  and  gear  box  driven  PMG  and  superconducting  generator.  In  this 
case,  four  PMG's  are  used  to  supply  the  required  electrical  power. 

Before  comparing  system  point  designs,  it  is  useful  to  point  out  a 
few  items  of  interest  in  Table  5.  In  general,  the  PMG  is  a  high  speed, 
high  frequency,  and  low  output  voltage  device;  the  superconducting  gen¬ 
erator  is  just  the  opposite,  i.e.,  low  speed,  low  frequency,  and  high 
output  voltage.  These  features  should  be  carefully  considered  when 
selecting  a  machine  for  an  application  where  these  parameters  are  critical. 
The  PMG  is  much  smaller  in  volume  and  somewhat  higher  in  efficiency  than 
the  superconducting  generator.  The  superconducting  generator  is  much 
lighter  than  the  combined  PMG's  at  the  higher  power  level;  however,  at 
the  lower  power  level,  the  combined  PMG's  are  lighter  if  the  rotational 
speed  is  high  enough. 

A  summary  of  the  system  point  designs  is  given  in  Table  6.  No 
volumes  are  given;  however,  generator  volumes  are  given  in  Table  5. 

System  weight  relative  to  the  lightest  system  for  each  of  the  two 
power  levels  is  shown  in  the  right  most  column.  Note  that  the  gear  box 
driven  PMG  is  the  lightest  system  at  the  10  MW  power  level  and  the  directly 
driven  superconducting  generator  is  the  lightest  system  at  the  20  MW 
power  level . 
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Plots  of  component  weights  in  Figures  25  and  26  enable  one  to  more 
easily  see  how  the  weight  distribution  within  the  system  changes  as  the 
generator  design  changes.  Note  that  the  fuel  weight  dominates  at  both 
power  levels.  A  14,900  RPM  turbine  using  JP4-L0X  fuel  with  a  specific 
propellant  consumption  (SPC)  of  7.0  lbs/(horsepower-hr)  is  used  for  the 
10  MW  power  level,  and  a  10,500  RPM  turbine  using  JP4-L0X  fuel  with  an 
SPC  of  7.0  1 bs/(horsepower-hr)  is  used  for  the  20  MW  power  level. 
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Figure  25.  System  Weight  Distributions  for  the  10  MW  Point  Designs 
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POINT  DESIGN  NO.  5  (PMG) 


POINT  DESIGN  NO.  6  (SUPERCONDUCTING) 


POINT  DESIGN  NO.  7  (SUPERCONDUCTING) 


Figure  26.  System  Weight  Distributions  for  the  20  MW  Point  Designs 
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SECTION  VI 
CONCLUSIONS 

Equations  for  calculating  required  generator  and  turbine  powers  were 
derived  assuming  that  the  component  efficiencies  in  the  system  were  known. 
Simplified  weight  and  volume  algorithms  for  high  power  permanent  magnet 
and  superconducting,  synchronous,  electrical  generators  were  developed 
and  used  to  obtain  parametric  results.  These  algorithms  are  in  the  form 
of  computer  codes  which  feature  simple  input  format  for  the  important 
input  parameters.  This  allows  one  to  make  use  of  the  algorithms  for  other 
parameter  ranges  of  interest  to  obtain  weight  and  volume  estimates.  In 
addition,  generator  efficiency,  which  is  a  very  important  system  parameter, 
is  calculated.  Further,  it  was  found  that  fair  comparisons  of  generator 
designs  (of  the  same  class)  at  different  power,  RPM,  or  voltage  ratings 
can  be  made  only  if  the  specific  electric  loadings  for  the  different 
designs  are  identical. 

General  conclusions  regarding  the  two  kinds  of  electrical  generators 
include:  The  PMG  is  a  high  speed,  high  frequency,  low  voltage,  low  volume 
machine,  while  the  conductively  shielded  superconducting  generator  is  a 
low  speed,  low  frequency,  hiqh  voltage,  large  volume  machine.  The  weights 
of  the  two  generators  are  not  very  different  at  10  MW;  however,  the  super¬ 
conducting  generator  is  significantly  lighter  at  the  20  MW  power  level. 

The  PMG  has  a  higher  efficiency  than  the  image-shielded  supercond'icti ng 
generator  at  power  levels  below  50  MW.  Volume  estimates  for  image-shielded 
superconducting  generators  are  more  than  three  times  greater  than  for  the 
PMG;  however,  it  must  be  noted  that  the  stator  coolant  is  integrally  built 
into  the  frame  of  the  superconducting  generator,  while  the  coolant  for  the 
PMG  requires  external  storage  and  heat  exchanger  space.  This  offsets  the 
volume  advantage  of  the  PMG  to  some  degree,  although  the  exact  volume  of 
the  external  coolant  loop  for  the  PMG  has  not  been  calculated.  Also, 
image-shielded  superconducting  generators  are  inherently  of  large  volume 
but  very  low  weight.  Superconducting  generators  using  iron  shields  will 
have  much  smaller  envelope  volumes;  however,  the  iron  will  tend  to  make 
these  machines  heavier  than  their  image  (conductively)  shielded  counter¬ 
parts.  Although  these  general  conclusions  are  valid,  a  true  comparison 
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of  the  two  classes  of  electrical  generators  can  only  be  made  by  comparing 
complete  power  supply  systems  as  described  below. 

The  generator  algorithms  plus  the  component  algorithms  in  Appendix  C 
were  used  to  obtain  system  point  designs  at  10  and  20  MW  power  levels. 

From  the  system  point  designs,  the  two  classes  of  electrical  generators 
can  be  compared.  From  Figure  27,  it  can  be  seen  that  the  PMG  results 
in  the  lightest  system  at  the  lower  power  level,  while  the  superconducting 
generator  results  in  the  lightest  system  at  the  higher  power  level.  The 
two  lightest  systems  are  quite  different,  besides  having  different  kinds 
of  generators.  The  PMG  is  gear  box  driven  and  operates  at  a  high  electrical 
frequency,  but  the  superconducting  generator  is  directly  driven  from  the 
turbine  and  operates  at  a  lower  electrical  frequency.  These  differences 
account,  in  part,  for  the  relatively  small  variance  in  system  weight 
between  the  point  designs  at  a  given  power  level.  For  example,  even 
though  the  superconducting  generator  is  much  lighter  than  the  four  PMG 
modules  at  the  20  MW  power  level,  the  low  electrical  frequency  of  the 
superconducting  machine  imposes  a  weight  penalty  in  the  power  conditioning 
subsystem  (transformer  and  filter).  The  weight  differences  between  the 
lightest  and  second  lightest  systems  are  440  and  more  than  1500  pounds 
for  the  10  and  20  MW  systems,  respectively,  and  are  significant  from  a 
mission  payload/performance  perspective.  Although  the  individual  compo¬ 
nent  algorithms  are  only  10%  accurate,  a  comparison  of  systems  using  the 
same  algorithms  should  be  valid.  In  other  words,  even  with  the  most 
detailed  point  designs  for  every  single  component  of  the  power  supply 
system  (rather  than  the  approximate  algorithms  of  this  report),  the 
relative  weight  standings  shown  in  Figure  27  for  the  seven  point  designs 
should  not  change  dramatically;  however,  the  absolute  magnitudes  of 
the  system  weights  would  most  likely  shift  upward. 

The  lower  efficiency  of  the  superconducting  generator  is  offset  to 
some  degree  by  the  absence  of  a  high-power  transformer  and  the  attendant 
gain  in  system  efficiency.  This  is  illustrated  by  the  *=>ct  that  the  fuel 
supply  for  both  PMG  and  superconducting  systems  weigh  about  the  same. 
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This,  and  the  fact  that  the  weight  of  the  fuel  is  a  large  part  of  the 
total  system  weight,  are  other  reasons  for  the  relatively  small  variance 
between  the  weights  of  the  system  point  designs  at  similar  power  levels. 
However,  again,  the  importance  and  significance  of  the  weight  differences 
must  be  emphasized.  Also,  the  higher  operating  frequency  of  the  PMG  tends 
to  make  the  power  conditioning  equipment  lighter  for  this  system;  this  is 
not  apparent  in  Figures  26  and  27  because  the  PMG  power  conditioning 
subsystem  includes  a  high-power  transformer  and  the  superconducting  power 
conditioning  subsystem  does  not. 

There  are  certainly  factors  other  than  system  weight  which  must  be 
considered  in  selecting  a  generator.  For  example,  in  point  design  number  6, 
which  is  the  lightest  system  at  the  20  MW  power  level,  the  superconducting 
generator  is  directly  driven  by  the  turbine  at  10,500  RPM.  Th’s  will 
require  careful  design  of  the  superconducting  field  coils,  coil  support 
structure,  and  helium  management  system  for  operation  at  this  high  speed. 

A  12,000  RPM  superconducting  rotor  has  already  been  built  and  tested  by 
Westinghouse  Electric  Corporation  (Reference  14). 

It  would  obviously  be  desirable  to  improve  the  accuracy  and/or 
flexibility  of  the  two  generator  computer  programs.  One  area  for  improve¬ 
ment  would  be  to  more  accurately  determine  the  weight  and  volume  penalties 
associated  with  interconnecting  the  5  MW  PMG  modules.  Another  would  be 
to  expand  subroutine  MFLD  in  the  superconducting  generator  program  to 
include  other  field  winding  configurations.  In  addition,  it  would  be 
useful  to  consider  the  effect  of  reducing  the  "air  gap"  in  the  super¬ 
conducting  machine  for  higher  rotational  speeds  as  the  rotor  decreases  in 
diameter.  This  could  be  done  by  reducing  the  thicknesses  of  the  EM  shield 
and/or  the  torque  tube.  Calculation  of  the  machine  inductances  for  dif¬ 
ferent  designs  in  order  to  determine  voltage  regulation  and  the  effects  of 
operating  with  rectified  output  or  pulse  forming  subsystems  is  also  needed. 

From  the  overall  system  viewpoint,  estimates  of  component  volumes 
such  as  the  turbine,  fuel  supply,  cooling  subsystems,  and  power  condition¬ 
ing  components  are  required.  One  further  step  would  be  to  consider 
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packaging  and  interconnection  constraints  in  assembling  the  power  supply 
system  into  the  compact  quarters  of  an  aircraft. 

Computer  programs  such  as  those  described  in  this  report  are  used  as 
part  of  the  AFAPL  computer-aided  power  system  design  in-house  research 
effort.  These  programs  are  used  to  predict  the  weight  and  size  of 
components  in  high-power,  airborne,  power  supply  systems.  These  programs 
yield  approximate  weights  and  volumes,  since  they  are  not  detailed  design 
programs;  the  results  should  nevertheless  be  useful  in  feasibility  and 
tradeoff  studies  concerning  these  two  kinds  of  generators.  One  of  the 
major  advantages  of  this  kind  of  "approximate"  program  is  the  fact  that 
it  adapts  easily  to  "computerization"  and  allows  a  large  number  of  point 
designs  to  be  made  with  moderate  computer  time. 
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APPENDIX  A 

GENERATED  VOLTAGE  IN  ALTERNATING-CURRENT  SYNCHRONOUS  MACHINES 

Figure  28  illustrates  the  rotor  of  an  AC,  synchronous  generator 
(alternator)  with  the  projection  of  stator  coil  (aa')  area  shaded.  All 
dimensions  are  assumed  to  be  in  meters  and  the  angle  0  is  measured 
as  indicated.  The  area  of  the  projection  is  ttDL/p  square  meters  for  a 
p  pole  machine.  It  is  assumed  that  the  rotor  produces  a  sine  wave  of 
magnetic  flux  density,  which  is  shown  in  Figure  29  at  a  given  time,  t. 

The  average  flux  density  linking  coil  aa'  at  this  time  is  given  by 

/TT-Ult 

B  sine  d6  =  -^  cosut,  (webers/m2)  (76) 

-U)t 


where  0  is  the  angular  displacement  and  B  is  the  peak  value  of  the  sine 
wave  of  magnetic  flux  density.  The  average  flux  linking  coil  aa'  at 
time  t  is  just  the  area  times  the  flux  density,  or 

<l>avg(t)  =  ‘f’gCOSwt,  (webers)  (77) 

where 

<J>e  =  2DLB/p  (webers)  (78) 

is  the  total  "effective  webers  per  pole"  and  can  be  written  as 

4>e  =  [ttDL/p]  x  [2B/tt] ,  (webers)  (79) 


where  the  first  term  is  the  area  projected  onto  the  rotor,  and  the  second 
term  is  the  average  value  of  a  positive  half-cycle  of  a  sine  wave. 

For  a  coil  with  N  turns  connected  in  series  and  linking  the  same 
flux,  the  generated  voltage  is  given  by  (Reference  10) 

Eg  ’  '  N  Ht  *a«g(t)  (voUs)  (80) 

or 

Eg  =  -  N  (d<t>e/dt)  coswt  +  iuN<}>e  sintot  (volts)  (81) 
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If  <j>  does  not  change  with  time  (as  is  usually  the  case). 

Eg  =  uN<|)e  si  nut.  (volts)  (82) 

The  rms  voltage  per  phase  is  just 

Vph  =  ojN  y/T 

Vph  =  ^  N<J>e  (volts  rms)  (83) 

By  including  a  winding  factor,  Kw>  which  accounts  for  the  distribution  of 
the  stator  coils  relative  to  the  spacing  of  the  field  poles,  the  generated 
phase  voltage  becomes 

Vph  =  4.44  f  1^  N  [2DLB/p] ,  (volts  rms)  (84) 

where  the  expression  for  N  is  given  below  Equation  9  in  the  text.  Con¬ 
verting  from  MKS  to  English  units  then  yields  Equation  9  in  the  text. 
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APPENDIX  B 

COMPUTER  CODE  LISTINGS  AND  SAMPLE  OUTPUTS 

The  following  computer  code  listings  and  sample  outputs  are  given 
for  reference.  Sample  outputs  were  obtained  directly  from  the  listings 
indicated.  Fortran  IV  programming  language  was  used  throughout  and 
comment  cards  were  liberally  employed  which  should  make  the  programs 
easier  to  understand. 

Input  parameters  appear  at  the  beginning  of  each  program  and  can 
be  varied  by  changing  the  value  of  any  input  parameter  at  this  point 
in  the  program. 

All  of  the  component  mass  densities  for  the  PMG  and  superconducting 
generator  are  given  in  Tables  7  and  8  respectively. 

Figure  30  is  referenced  in  the  superconducting  generator  program 
SCGEN. 


TABLE  7 

DENSITIES  OF  PMG  COMPONENTS 


COMPONENT 

DENSITY  (lbs/ 

Stator  yoke 

0.2815 

Stator  teeth 

0.2815 

Stator  insulation 

0.100 

Rotor  (overal 1 ) 

0.238 

Copper  conductors 

0.320 
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DENSITIES  01  SUPERCONDUCTING  GENERATOR  COMPONENTS 


COMPONENT 


Environmental  shield  (components 
1  and  2  in  the  computer  program) 

End  bells 

Field  coils  (epoxy-impregnated  modules 
with  a  75%  copper  packing  factor) 

Anti-drive  end  structure 

Drive  end  structure 

Torque  tube 

Bore  tube 

Field  winding  support  structure 
Electromagnetic  shield  (rotor  shield) 
Copper  conductors 
Stator  insulation 
Bore  seal 


DENSITY  (Ibs/cu.in. ' 


0.75X  0.32+  0.25X 
0.10  =  0.265 
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PERMANENT  MAGNET  GENERATOR  WEIGHT  AND  VOLUME  COMPUTER  PROGRAM 


°ROGR A  1  ^MGI TUI»  IT ) 

)l PENSION  ->P<510) 

C  This  PROGRAM  C Al  CUL  4  TtS  THE  IEIGhT  An  /OLUN  *  Op  p-RMAME  <T  MAGN_T  GE  NtMTORS 

^  **«*.*..*  .........  **..♦.*..  .........  •***«*..« 

c  ...  PERMANENT  MAGNET  G.H"!-TA»  OAST  LI  I  DESIGN  *  *  * 

*****  *****  *****  *****  *****  *****  *****  ***** 

*****  TH-  r  ILL  0  4 1  IG  GENERATOR  HIRAM.T.RS  MOST  ’>;  Sp_CIFj.-0  TOR  ♦♦*»**»* 

A  Gx  IE RAT  OR  POINT  CESIG  I 
>PrCTFY  TrC  GENERATOR  ONTO  JT  POWER  IN  WATTS 
JG-S. 0  3 "*Gj 

C  GENERATE!  *  RP  1 

RPM= Id  J  0  3 • 

C  PHAS'  VOL  TAG-  (VOLTG  RIG) 

vph=s« a.  ■ 

C  LOAD  PONx  R  r  APTO 

PR=.32 

C  COPPER  C  JRRr  NT  1”  NS  1 1  Y  II  TH.  STATOR  ( A  1»S /S  J.  Ill .  ) 

OJ=3G273,3C 

C  ROTOR  II  -»  SPEED  (FI/ SEC) 

VT -h 23  «  3  E- 

C  NUTTER  OF  POlES 

P=1A. 0 

C  NUHR„P  OF  phases  in  th;  stator  WINDING 
UH  =  1. 

C  *•*  DEFINE  HE  PITCH  i]  A  FRACTIJI  OF  PI  (OR  IRC  D-GRE.S  ELECTRICAL) 

RI=3. 1P1 SORT 
PTTCH=(R. 1/6.0) *°I 

0  .**  SPECIFY  TH  :  STATOR  INSJlATIJN  RATING  IN  VOlTS/'-IL 
vP'iiL=nj.ct 


C  PEAK  MAGNETIC  FL'IX  DENSITY  TN  TH-  AIR  uA»  (LIN..S/SO*  IH.  )  UlOER  LOAF 

TG  =  *  ’**00.  3 

C  SPACE  F&CTOP  FOR  STA'OR  TAR  COOLI  Ij 

F  =  *  8  3 

C  RATIO  OF  rUU  J’NSITY  IN  THE  TCOTH  TO  MaT  IN  THE  AIR  GAP 

ALPH=2.2J 

C  -<ATIO  QP  SHAFT  0)  TO  ROTOR  V) 

FR- 0**06 

c  :no  tu=i  angle 

CH I =P 1/ . *  6 

c  *»*  sr t  the  allowail:  fljy  o  :nsity  in  thl  rack-iron  (ltnes/in*»2)  ***»♦♦♦*»♦• 

■»YO«E  =  10303C.G 

C  *.***.*....**..*  i;Nj  jF  D)T4T  OESIG j  specification  ................. 

C  ***  MAXIMUM  SPPCIFIC  ELECTRIC  LOADING  (A“PERC-CONDUCTORS/INCh) 

DMAXsAOIO. 

C  RATIO  OF  FLU <  DENSITY  IN  THE  GA°  TJ  THAT  IN  THE  STATOR  YOKE 


HGGY-DG/  HY  OX  E 

C  ***  CAoCULATE  The  P-AK  FlUR  JEN5ITY  IN  THx  YOKx  IN  KLllxS/IN*»2  **« 
TYKL  =  FJYO<  E/1  00  j.  C 

C  •••  CALCULAT.  The  RPAK  F.UX  DENSITY  IN  T Hi  TEETH  IN  Kt IN ES/ IH**<?  ♦** 


)T'th-tg*alph/io  ;  c. o j 

C  *«.*****»*.  SlT  T ME  COMPONENT  MASS  DENSITIES  (IN  LUS/  CJDIC  INCH) 
C  JFNSITY  OF  STATOR  YO<E 

RY=.231S 

C  density  of  stator  t.eth 
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RT  =  . 2815 

C  DFNSITY  OF  INSULATION  IN  T Hi  STATOR  SLOTS 

RINS=.19 

C  ROTOR  OFNSITY 

RR=.23S 

C  COPPER  DENSITY 

RCU  =  . 32 

C  »*»  BEGINNING  OF  GENERATOR  DESIGN  LOOP  ********** 

C  INITIALIZE  THl  NO.  OF  CONOUCTORS  PER  SLOT  TO  TWO 

CS  =  2. 

C  ANS  =  NO.  OF  SLOTS  HER  POLE  PER  PHASE 

ANS=4. 00 
•550  CONTINUE 

C  CALCULATE  THE  OUTPUT  POWr R  IN  MEGAWATTS 

W=1.S-06»PG 

C  *♦*  oALCULATE  THE  FREQUENCY  (HoRTZ)  ********** 

FREQ=(P/2.C>MRPM/oC.00> 

PRINT  ITT, VPH,RPH,W, FREQ 
C  ******  CALCULATE  THE  ROTOR  3IAHETER  ******** 

Q=(T20.*i/T)/(PI*RPM) 

C  ***  CALCULATE  THE  QUANTITY  L  PRIME  USED  IN  THE  COMPONENT  VOLUME  EXPRESSIONS 
ALP=(0/ <1.h14*SIN(CHI ) t ) *SQRT (1. -COS  <2.*PI/P> ) 

C  CALCULATE  THE  THICKN'SS  OF  THE  STATOR  YOKE ( RACK  IRON) 

)Y  =  ( D/P >  4  OGOY 

C  CALCULATE  THE  RATED  PHASE  CURRENT (AMPS  RMS) 

AMP=  <PG/PF)/(PH*VPH) 

C  01  = INSULATION  THICKNESS  Fo R  LINE  TC  NEUTRAL  VOLTAGE  (WALL  INSULATION) 

J1  =  ( SORT ( 2. ) » VPH/VPMID  *  1. E-u  3 
C  J2=INSULATI0N  THICKNrSS  BETWEEN  COILS 

02=01*2. j 

C  CALCULATE  THE  TOTAL  to.  OF  ACTIVE  oONOJCTOkS  IN  THE  STATOR 

664  7=ANS*P*PH*C3 

C  CALCULATE  THE  TOTAL  10.  OF  SLOTS 

5L0T3=ANS4P»OH 

C  **********  CALCULATE  TH  :  P0SSI5LE  NUMBER  DF  PARALLEL  path?  ************* 
c  ********  JN  THE  STATOR  WINDING  ************** 


453 


L5C 


45? 


452 

45 

455 
457 

********************  pjoalLEL  PATH  CALCULATION  ********************** 

»•*  CALCULATE  THE  WTNUINS  FACTOR  As  THc.  PRODJCT  OF  THE  PITCH  FACTOR 


CTG  =  A  NS  *  3 
DIV=1. 0 
1  =  1 

DIVR=FIG/DIV 
PDT  =  IFIX  (  0  IV  R) 

PPT2=0IVR-PPT 
IF (PPT2. EQ.C.O)  GO  TO  453 
GO  TO  455 
°P ( I ) =OIVR 

:h=i 

T  =  I » 1 

1IV=0IV* 1 . 0 

IF (OXV.GT.FIG)  GO  TO  452 
GO  TO  453 
CONTINUE 
PRINT  453 

CQRM A  ?  <1H1,1-»HPA\ALL“L  OATHS,/) 
DO  *»56  1  =  1, It 
OPINT  457,1,00(1) 

FDPMAT ( I1D,F2C . 7) 


-V:  ;  C 


r-°DDlc 


CBs>?*or  / 
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C  AND  THE  3  1ST ri jo t i on  factor 
c  ♦  calculate  th:  pitch  factor 
AKP=$IN(PITCH/2.00> 

C  *•*  CALCULATE  THE  DISTRIBUTION  FACTOR  FROM  THE  NUMQFP  OF  COILS  IN 
C  ***  .ACH  PHASE  GROUP  A  NO  the  Elcct-ICAL  PAJIANS  DcTWEFN  SLOTS  (  ^  ) 
3-PITCH/ (ANS»PU) 

AKU=( 1. 02/ANS) *  tSlNl ( ANS*D) /2. Cn) ) /SIN (9/2. 03) 

C  ***  CALCULATE  the  COMBINED  WINDING  FACTOC 

akw=ak°*  ako 

C  »**  ITFRATi  ON  PPATH  FOR  GIV  E N  VALUES  OF  AN3  ANO  CS 
00  bo  8  1  =  1,11 
PPATH=PO( I) 

C  ***  CALCULATE  the  CUsR: NT  °E  R  CONDUCTOR  (ARPS  RMS)  »»»* 

AMPC-=AMP/PRATH 

C  ***  CALCULATE  THE  STACK  LENGTH  *♦* 

C  ***  THE  FACTOR  (CS/2. 00’) ** (-, 36S46)  ACCOUNTS  FOR  SLOT  LEAK AGt  AS  THE  NO. 
C  »**  OF  CONDUCTORS  P"R  slot  increases 

AL 1 VPH/  (  /  .  -.C  183-10  ♦°_>M»A<W»  (  (  P/PPA T H )  ♦  AN S*  ( CS/ 2 .  ) ) *0*Db*  (CS/2. 00  u  > 

&**  (-.  06S-.G)  ) 

C  *♦*  CALCULATE  the  specific  electric  loading  *»* 

QL= (Z*A1PC>/ (“I* (0  +  . S?>  > 

IF(DL.GT.T-IAR)  go  to  1?Z 

C  CHECK  PRASr  VOlTAG~(\/PH1  SHOULD  E QUA u  /PH) 

VPH1=3. *  339L-13* AKW*RP 1* J*AL»3G*<(CS/2.3D0>**(-. 0bS4o) )  *<  12/ ( PPATh 

L*PH) ) ) 

G  CHECK  POWER  <WC  SHOULD  "IUAl  W) 

WC=  AMP* >f*ph*V°H1 

C  CALCULATE  THE  SLOT  WIDTH (Ini 

S=(PI*0»( ALPH-1. ) ) / (ALPH* SLUTS) 

C  C=CONDUCTDR  (OAR)  Wl)TH  (INCHES) 

C=S-2.»01 

C  CHr C K  TO  SEE  IF  Th^  CONDUCTOR  WIDTH  IS  POSITIVE 

IF(C.LE.I.C)  GO  To  2jC 

c  **♦  calcjlat.  the  slot  depth  *♦* 

jT=(CS*A1PC>/(‘r*C*CJ)*(0;>-2.ij)*.H*Jl*2..*02 
C  »**  CALCJLAT-  THE  A  PE  A  Oc  A  SINGLE  CONDUCT  DR  IN  SQUARE  INCHES  *** 
AFrAMPC'C J 

WEIGHT  CALCULATIONS  (  EL  ECTRU'IAGNET  IC  WEIGHT  ) 

NY  =  DA  CK  IRON  (YOKE)  WEIGHT  (POUNDS) 

NY=RY* <°I*OY* ( J  +  2. *  IT ♦ OY ) * Al) 

C  WT  =  WEIGHT  OF  Trtr  T-ETH  t POUNDS ) 

NTrOT*(PI*(0*0T>* 07 *AL)/AL PH 
C  WTNS=  I (SULATIDN  W.ISHT  IN  THE  ST  A T OP ( P JJN CS ) 

WINS=RINS*(S*DT- (CS*AMPC> /SJ)»SL0T3*(AL*2.*ALP) 

C  WP=  ROTOR  WEIGHT (POUNDS) 

Wr-RR*( (JI/4. ) * ( 5**2) *( 1 • -FR*»2> ) *AL 
C  WCUrWEIGHT  of  T  H„  COPPER  IN  THE  STATOR  C 01 L S ( F OUNDS ) 

WCU=RCU* jLJTS»aC»CS* ( AL+ALP) 

C  WTOT  =  TDTAL  f LECT ROMA  GN. TIC  WEIGHT  (POUUOS) 

WTCT  =  WY*-WT*WIN3«-WR*W0U 
C  TOTAL  MACHINE  WIGHT  =  W  'ACH(PCUNJS) 

WMACH=1.20C*WTJT 

c  calcjlat:  the  s°:ciFrc  w  ight  in  lds/kw 
sw  =  (wmach/pg)  *i‘,jc.j: 

C  CALCULATE  the  frame  RADIUS ( IN) 

(FRM= (J*.5a*2.*OT*2.*nY*2.jG)/2. 8j 
C  FDD  I A  ?F  »  VI '  DH-*ET.R<I  )) 

FPD I A  =  2 • S  *  °FRU 


SsSSSSP 

erspkotwcs  legibly  . 
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C  'AlCULAT:  THa  total  IAJHINa  VOLUMa  If.  CJ-il  '  F.'T 
VOL  = ( PI* RFkM*»2> *  UL*Ald*OOS(CHI>  > /1/23, 3C 
c  »»  m-  a \/g .  t£mp-:ratjp.  >f  t,<-  stator  sajs  ia  assum.j  to  ol  tc  j •  g p : :  3  c: 
Tr=232.9 

C  »"  CALCJLAT.  Th£  STATOR  OAR  R-.  i  I  $  T  I V  T  T  Y  Ii  OH 1- INCHES  *♦* 

PHOC=o.  77?a-J7*(23*.5*T3>  /25l. 5 

C  ***  CALCULAT  £  T  Ha  OHHIC  HaATINa  Pak  CLHIC  INCH  OF  STmTO°  COPPFR  **♦ 
°IRV=Cj*aJ*RHOC 
RRL  Oa  =  P I R  V*  HC  J/RCU) 

c  *•*  talculat;  r h ..  hyst.r-sis  loss  in  th:  Yo<a  •  ** 

PHY  =  7.1£-95*FR:Q*0Y<L»*1.6,. 

C  *»*  CALCULATa  THa  £93Y  CJRRINT  LOSalS  In  TH£  YC,Ka  *♦♦ 

'A£Y  =  4.3r-09*fPr.  3»*2*3Y<L»*2 

C  »**  CALCOLAT.  TH.  HYSTSR-SIS  LJSSaS  IN  THa  T.aTn  »*♦ 

PHT  =  7.1£-0F*FP-:3*«T£Th**l.o0 

c  ♦♦♦  cucolat:  th£  ;qoy  current  lossss  in  th:  u.th  •** 

PFT=4. 3a-0  3*FR£0**2*  3TSTH**  2 

C  ***  CALCOLAT.  THa  LOSSES  10  WATTS  9Y  MJLTIRlYInG  3Y  TH.  WFIGhTS  IN  POUlOa 
r>YK  =  KY»  (  PHY  +  PSY  ) 

PTH=  W  T • ( PHT*P£T ) 

0  *»*  JAlOULAT-;  TH.  TOTAL  IRON  LOSSaS  *»» 

PIRON=PY<*PT  i 

C  *•*  aALCOLATI  THa  TOTAL  -TAT  3R  LOSS  (WATTS ) 
pLOS3=PRlOS+pIRON 

C  ***  aALCULATJ  th:  NACHIN.  aFFICItNCY  3  A  a  a  0  3  4  STaTGR  10SS£S  GNLY  *** 
aFF= (PG/ ( PG*  3L OSS) >*101.30 
GO  TO  665 

150  PRINT  151 

PRINT  1  >  a,  PPAT-I,  ANS,  CS,  TL 
GO  TO  663 
26 C  PRINT  251 

PP  TN T  16a, PPATH , ans,cs, ol 
.0  TO  aa  3 

<>6r  PPINI  30, SW,  RPH,  AMPC,  VPH1 ,  A  IP,  AKW 

PRINT  203, ANS.CS.Pf AT H, SlOTS, S,OT,nY 
PRINT  30  0 ,WY,  WT , WINS, W  R, WCU, H TOT ,WNACH 
PRINT  400,7, 3,  AL  ,  VOL  ,  QL 
URINT  431,  '  FF,RLOSS,PRLOS,RY<,PTH 
*6“  COnTINUT 
36  0  CONTI N','£ 

34  rORMA  T (//,  3X*CJ"C«  W 3 ISHT * 5X* RF  l*9X * AMPC *9 X*C ALCdL AT: J  PHASS  VOLTA 
<.„a*2X*CJ  (RaN T  (  AM°S >  »  2  X* 'W I  N  >NG  PTCH*OIST  F A CT  OR*  /  3F  1 2 . 3 ,  F  24 .  b  ,  2F2  .  . 
L6/1 

2C0  FOR  HAT (*  NO.SL/P/PH  03  PPATHS  NOM3SR  OF  SLOTS  SL 

LOT  WI0TH  SLOT  3£°Th  YOKF  Oa FT H* /3 F 1 0 . 3 , 5 X , *F 1 5 . 3 / > 

3C0  FORMAT (3 X*Wl IGHTS-L 3$ . * 4 «.» f  OKa* 1GX* TFETH*5 X* INSOL AT  I  ON* 7 X* KOTOR*  3X 

L  *COPP  rR* 3X*T JTAL*9X*HACHIn:*/1.X,7F15.2/> 

400  FORMAT ( 3X* TOTAL  NO.  OONO.  ROTOR  3IAUNI  STATOR  Lat(IN)  VOLUMl 
L*6X*SP£C.LOA0»/5Fl?.*/> 

401  FORMAT (/,5X,13HrFFICItNCY  =  , F 7. 2 , 3 X , ?HL OSS  =  ,f a . 1 , 3X , 3H0HM IC  =  , 

L  FM  • 1, 3X , 7HYC<F  =  ,F», 1 , 3X , AHTEiTH  =  ,F9.1,/> 

151  FORMATflorl  lOA)ING  £XCr.0£0) 

251  FORMA  T (1 3H  NO  ROOM  FOR  COPPLR) 

16b  F0°NAT(aF15.4,/> 

177  FORMAT (1H1 ,5H  V°H= ,FM. 2, 2X , *HRPM= , F 1 3 . 2 , 2X , 6HP0W £R  = , FI 2 . 2 , 2 X , 5hF R £ 
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LP*1-  HOC,..} 


5.00  FRLu*  21CQ.0* 


P ARALLr  .  PAT  hS 


56. 3  J  0  a  J  J  3 
29. 9  "09  3  30 
14. 3CG.J-C 
8.  C0C30J3 
7.  3  <10  JO  1  3 
4.  3  JGO  ,0 

2. orocjoo 

1.9003030 


SPEC.  WIGHT  *PR 

.J75  I63cr.30? 


CALCULATED  Phase  VOLTAGE  CURRENT (AHPS)  winong  ptch^oist  factor 
569.3G0030  J45o. 667220  .9  37424 


nO.SL/P/ph  CS  PPAThS  NUMBER  OF  SLOTS  SLOT  WIOTH  SLOT  DEPTH  YOKE  11PTM 

%•  0  30  ?.0EC  56. 633  To5. 0  33  .  381  .10-  ,27G 


ME  I jHl S-L9S ■  YOKE 
?51.3 J 


INSUL  AT  ION 
1C. 9b 


TOTAL  NO.  con.  ROTOR  01 A  C I  4)  STATOR  LENdN)  volume  spec. load 

JS6.93CC  7.96.1  12%. 71,1  7.2445  779.9o92 


EFFICIENCY  -  95.33  LOSS  =  245171.7  OHMIC  =  119C28.J  YOKE  *  1C6921.5  TEETH  =  19221,9 


S-fcC.  WEIGHT  “P*l  AHfC 

.194  l40CG.33-3  123.462 


CALCULATED  PHASE  VOLTA&E  CURRENT  ( AMPS )  NINONS  PTCH*OIST  FACTOR 
589.GGJGJJ  3455.667220  .937,2. 


NO.SL/p/Ph  CS  PPAThS  NUMBER  OF  SLOTS  SLOT  NTOTH  SLOT  DEPTH  YC*P  DEPTH 

4.3i«  2.000  28.030  168.000  .C8l  .166  .270 


NEISHTS-L9S.  YJKE 
127. 36 


INSULATION 

7.57 


TOTai  NO.  COnJ.  ROTOR  JlAllN)  c ATDR  LEN(lN)  VOLUME  SPEC. LOAD 

316.9303  7.96*1  62.3570  3.7645  1559.9385 


EFFICIENCY  *  95.  T1  LOSS  3  191  330.4  OHMIC  =  121601.1  YOKE  =  5-2.5. 9  TcETh  =  15481.3 


SPEC.  WIGHT  opm  AhpC 

.1C-  19030.000  2-6.905 


CALCULATED  PHASE  VOLTAGE  SU*R-:nMAHPS>  N I  HD  N  S  PIDhOIST  FACTOR 
588.000010  3456,667223  . 937-2- 


NO.  SL/P.PH  CS  “PATHS  NUMBER  OF  SLOTS  SLOT  MlOTH  SuOT  DFPTh  YOKE  DEPTH 

4.033  2.  030  144-3C  168.000  .G8l  .293  .2'-? 


NtlSHTS-LRS.  YOK. 


INSULATION 

6.11 


MACHINE  I 
622.16 


TOTAL  NO.  CON  9.  ROTOR  O  I  A  4  I N  )  STATOR  l£N(IN)  VOLUME  SPEC. LOAD 

316.030G  7.9641  11.1755  2-0239  3119. 877C 


EFFICIENCY  =  36.  74  LOSS  3  168-04.5  OHHC  =  12o752.7  YOKE  *  27938.  1  TEETH  *  13743,7 


LO/JING  E*CEOE) 
9.3300 


LOADING  E  *C 1 9£9 
7. 0000 


LOADING  E«C*;OEl 
4.0C0C 


LOADING  PACLDET 
2.3300 


LOADING  f  * C " DK T 
1.3300 


2.3330  5459.79-7 


2.300G  6239.7539 


2.00TC  1C919.569- 


2.030S  21839.1318 


2.0000  -3678.2776 


**/r, 
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SUPERCONDUCTING  GENERATOR  WEIGHT  AND  VOLUME  COMPUTER  PROGRAM 


C 


C 

c 


PROGRA i  SGGEN (OUT PUT ) 

DOJRLE  PRECISION  PMU,  RM1  ,KM2,RSM,0X,  OY,  XJ  ,  Y  0  ,  A  I 

DOUBLE  PRECISION  dRl,3Fl>9R2,3F2,BRM,3FM,3R,BTM,BRR,BFF,3T 

DIMENSION  VOL ( 13) ,P( 1 3) 

RlAL  <W,JC,MUC,L,LL,LC,IPHASt.,LPRlME,LSM 

****  ********  ........  ********  ******** 


♦**  SUPERCONDUCTING  GENERATOR  BASE  LINE  DESIGN  ******* 

........  »*♦*.**.  *....***  ........  ******** 


**********  THE  FOLLOWING  PARAMETERS  1UST  9E  SPECIFIED 
**********  FOR  a  uENcRATOR  POINT  DESIGN 
♦  ♦♦  SP-.CIFY  THE  GENERATOR  OUTPUT  POWER  IN  WATTS 
POWER.  -2fl.0E*06 

***  SPECIFY  THE  GENERATOR  RECTIFIED  OUTPUT  VOLTAGE  IN  VOLTS 

vo:=<fjiOi). 

.-PACIFY  THE  ROTATIONAL  SPEED 
HP  1  =  6  J  3  0  • 

***  SP'CIFY  THE  NUM3LR  OF  SLUTS  P_R  POLE  PER  PHASE 


ANS  =2  3  • 

...  SP  CICY  THE  NUMBER  OF  CONDUCTORS  PER  SLOT 


OS =6. 90 

***  SPECIFY  THE  EQUIVALENT  POWER  FACTOR 


************* 


PF  =  • Mo  E 

*»•  ^P.CIFY  THE  NUMBER  OF  PHAScS 
PHASu=3.Q3 

*•*  SPECIFY  T  H  STATOR  CONDUCTOR  CURRENT  OENSITY  In  AMPS/SO.  IN.  *** 
J 0  -  ? 3 L  :  G.  DC 

...  SPECIFY  The  SUPERCONDUCTING  FIELJ  WINDING  COIL  MCOULE  OVERALL 
CURRENT  DENSITY  IN  AMPS/CM* *2 
ACM-L5.00. 70 

iP’CIFY  THE  TIP  SPi_EJ  OF  THE  FIELD  COILS  IN  FEET/SEC  *** 


V  T -471.  1 

...  -;p  d IF Y  THE  DISTANCE  FROM  THE  OUTER  RADIUS  OF  THE  STATOR  WINQINi  TJ 
•  * »  Th  .  INNER  RADIUS  JF  THr  ENVIRONMENTAL  SHIELD  IN  INCHt.S 


S=7  .0 

...  SP  CIFY  THE  SUPERCONDUCTOR  MINIMUM  3l.N0  RADIUS 
3=1.593 


***  NOTE  THAT  THt  PROGRAM  PRESENTLY  CALCULATES  The  HEIGHT  (A)  OF  A  FIELD  COIL 
***  MODULE  AS  A  FIXED  FRACTIUN  OF  THE  RADIUS  OF  THl  ROTOR.  THAT  IS,  A=.32b7*RR 
THIS  CAN  3E  MODIFIED  TO  ALLOW  3PECIFICATIO  I  OF  ANY  FIELD  COIL  MODULE 
***  H.IGHT  3 Y  SIMPLY  DELETING  THIS  CARD  IN  THE  PROGRAM  AND  DEFINING  A  AT  THIS 
»**  POINT  in  THt  PROGRAM.  FOR  EXAMPLE,  IF  A  WlRE  TO  3E  2.30  INCH_S,  PLhCE  A 
♦♦♦  CA  RD  WITH  A=2.30  AT  THIS  POINT  IN  THt  PROGRAM. 


...  sp-CIFY  THE  ARMATURE  WINDING  DISTRIBUTION  FACTOR 
KW- 1.3 J 

***  3° .C I F Y  THE  STATOR  BAR  PACKING  FACTOR  <  I.E.,  THl  RATIO  OF  COPPER  TO 
*•*  INSULATION  IN  A  SINGLE  RECTANGULAR  SLOT 
FTAR-.Lpj 

o'j:CIFY  THE  DISTANCE  FROM  THE  EM  SHIELD  OUTER  RADIUS  TO  THl  DOPE  S_AL 
...  j-n-.R  RADIUS  IN  INCHES.  THIS  IS  THE  PHYSICAL  AIR  GAP 
G=. 2  D  C 

***  SP:CIFY  THE  STATOR  3AR  END  TURN  ANGLE  IN  RADIANS 
X  —  . u39 

***  HE  IS  TH'_  HEIGHT  AltO  WIDTH  OF  A  SQUARE  ,  HDLLOW  CONDUCTOR  WITH  A  WmLL 
***  THICKNESS  OF  TW  METuRS .  HO  IS  ALSO  IN  METERS. 
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SIGNIFICANT  NUMBER  OF  PAGES  WHICH  DO  NOB? 
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C 


H0*3.5J£-J3 
TI-H3-2. *TW 

**********  £N3  or  POINT  DESIGN  SPECIFICATION 


*********************** 


*****  THE  NUMBER  OF  POLES  IS  L I M I T r 0  TO  FOUR  FOP  THIS  PROGRAM  ******* 
POL_3  =  -».  0  0 

***  JgL  T  A3  =  T  HE  JORE  TUBE  THICKNESS  IN  INCHES 
DE-TA3*.!>? 

***  oeltat=the  tgroue  t joe  thickness  in  inches 

JEl  T  A T  *0.320 

DELTASONE  electromagnetic  shield  THICKN.SS  in  inches 

DELTAS*. 7s 

***  OE  _  T  A  E  =  T  rl  E  ENOBELL  THICKNESS  IN  INCHlS 
DE.-TAE*.5G 
PI  =  3. 1  ■*13927 
AN3=PI  /° Oi.ES 

♦**  JO  IS  T  Hi.  PEPHEABILITY  ok  a  vacuum 
Uj*1.2r>6o37E-0B 

•»*  SIGMA  IS  THE  CONDUCTIVITY  OF  ALUMINUM  IN  MHO 3/ Mi. TER 
SIGMA*  3.7Ytt-C7 
i  3  j  COlTINUE 

***  CALCULATE  THE  FIELD  HINGING  OD  *** 

BR=(72j.*VT)/(PI*RPM) 

*•*  RR  Ig  THE  RADIUS  OF  THE  FIELD  WINDING  STRUCTURl  IN  INCHES 
RR=0R/2. 

***  CALCULATE  THl  BORE  SlAL  THICKNESS  IN  INCHES  ************ 

V0CK* VDC/1 1  0  C  .  <j 
AK*1.7  32 

ei*-..  i: 

l0=2.11 

*”  ddsel  is  the  tore  seal  Thickness  in  inches 

JHSlL*  1 1.  0a3  -G-+*  <EI/ED)  *  (AK*VOCK)  ♦♦l.tfl) 

***  gBaR  IS  THE  INSULATION  THICKNESS  BETWEEN  STATOR  bars  OF  THE  SAME  PHASE 


»  *  * 
*  *  * 


*  *  * 
*  *  * 


IN  INCHES 
OiAR*. 12 Jj 

JPrt  IS  THE  INSULATION  THICKNESS  BtTWcEN  STATOR  BARg  OF  DIFFERENT 
PHASE  VOLTAGlS  IN  INCHES 
OPM*. FJO 

jCON  is  the  stator  conductor  insulation  thickness  IN  INCHES. 

DCOM*.  .  1-J 

lALCULATE  TH;  GlNlRAT  JR  line-t  0- neutral  volt ag^  in  volts 
VPM A gE=V DC/ 2. 3 

CALCULATE  TH.  h.IGHT  OF  a  WINDING  MODULE  *** 

S£“  MOTl  UNJ: p  GENERATOR  INPUT  PARAMETlR  DEFINITION 
A= ; . 32 >7*RK 

CALCULAT.  TH!  WIDTH  OF  A  FIELD  WINDING  MODULE  *** 
f  =  .  9*  (  ’.R-A*C03  (ANu)  )*SIN  (  ANG)  -  B 
^AiCULATE  TH-  AR.A  OF  A  MODULl  IN  CM**2 
ARCM=A*T*2.gL*2.FH 

CA_rjLATE  TM'  total  current  per  module  side 


CU  <R*  ACM* ARCH 
CALCJLATl  the  CURRENT 
AI  =  C')RR/l+4.  Cd 
kM* RK*. .234 
A1-A*. 1294 
TM=T*. 12g+ 


PEP  WIRg  BASED  ON  14*  WIRtS 
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.H-'S  UR  I  l  RM**  2-  (.  :  2  34*  3  +  T“)  **2  ) 

jx=  <  r/12. >  * .  :2s., 

JY  =  (4/12. )  *.(.374 
X j - . J  ?  4*'1*0X/2,C 

y,j=ch-am»-oy/2.  o 

2  **♦  OA.CULAT.  THl  to  OF  THE  ST  AT  OX  *** 

0*2.*  (RR*3ELTAT*0iLTAS*G> 
x* j/2.  :  o 

2  ♦**  CALC JLAT _  THC  TOTAL  NUMBER  OF  1L  OT  j  IN  T HE  STATOk 
oLOTS=4N3*PHAS-.*POLES 
2  ♦♦*  cAlCJLATE  THE  SCOT  WIJTh  IN  1NCH.S 

SWUM*  (PI*0-PHAS.*PGL  ES’UHH-  (  AMS-l.  C  )  ’PHAir»PCL2  3»03AP)  /  slots 
2  **»  CA-CJLAT.  THc  PHASE  CURRENT  IN  AIP.RcS(KrtS) 

IPHA2*=:PJNlR/  (PF*PHASE*VPHAS2> 

2  ***  CAL  COL  A  T  _  THt  AR.A  OF  A  SINGLE  STATOk  CONDUCTOR  IN  SQUkRE  INCHES 
AC  = IPHASS/JC 
SCHK=SWIDTH-2,C*OCON 
[F  (SChK.L,.  J.E)  «0  TO  Ipu 

C  *♦*  JA  lo  _  UAL  TO  THE  THICKNESS  OF  TH.  90Rc  SEAL  PLUS  THE  THICKNESS  **»» 

2  *»*  Or  TH  “  5  T  A  T  0  WIND ING  *** 

JA=  (2j  ’  I  PHASE)  /  (SCHK»FMAR*  JO  )  *  ( CS -1 .  C  )  *  .  0 1  *  2 . 2  »DC0N+ J  33  it 

c  *‘0  liA  t“  1 1  the  thickn.  ss  of  the  armatjr.-.  **» 

•;a  ri-OA-aHSFL 

r  ARM*  (  C  o* AC* SLOTS)  /(RI*J*DARM> 

c  ♦**  CAlCJL-T.  THl  SPECIFIC  .LsGTRIC  LOADING  (A  IP-2 ONOUCTOKS/I NCh) 

3=  ) A  RM  *  F  ARM*  JC 

2  ***  -IS  IS  THE  INSIDE  RAOUs  OF  The  ENVIRONMENT  AL  ihage  shield  *** 

P3  =  R-.*  licT  AT*  JELTAS  +  GV3A  +  S 

2  **♦  CALCULATE  Th:  INSIDE  RADIUS  OF  THC  STATOR  WINDING 
WHO  1- 2Rf  JCLTAT+OELTASfG  +  OPSEL 
2  »**  CAlCOL A  Ts  the  outer  radius  OF  THE  sTATOR  winding 
RH  32  =  PHO 1  *  J4-DC0N-03SEL 

2  *♦*  CONVERT  TH..  STATOR  WINDING  INNER  RADIUS  TO  M.T  iR? 

RMl  =  ."25-t*RK0i 

2  ***  CONVERT  TH,  STATOR  WIUING  OUTsR  RADIUS  TO  ML T EPS 
<M2=.  j  ?S7»PH02 

2  ***  PH  1  IS  T  HE  ANGLE  F  SON  TH-  COIL  DOPE  (  SEE  FIGOR.  30  IN  THIS  APPENDIX  ) 
P  10  =  3.  0 

2  *•*  CALCULATE  THE  lAOIAL <3k1 >  AND  AZIMJTHAl (DF1)  MAGNETIC  FLUX  DENSITIES  AT 
2  **♦  ThE  INNER  RADIUS  OF  THE  STATOR  WIN3ING 

C4.L  Mr LJ (RM1,PMU,PR1, OF  1 ,XC,Y  0,OX,OY  ,  A I ) 

2  ***  CALCULATE  THE  RADIAL  <3R1 )  AND  A Z I  HUT H A L  10F 1 )  MAGNETIC  FLUX  DENSITIES  AT 
C  •**  TH.  OUTER  RADIUS  OF  THE  STATOR  WINDING 

CALL  1FL0(RM2,PMU,BR2, 3F2,X0,Y0,DX,DY,AI) 

PRINT  s67,  3F 1 i OF  2 
7  K  7  FORMAT <//,2r23. 5,//) 
iRPl=DPl 
dR»2=P  <2 

C  ***  CALCULAT-  THE  AVERAGE  OF  THE  TWO  RADIAL  FLUX  DENSITIES  AND  CONVERT 
2  »**  TO  KIL  (GAUSS 

IK  >  =  A  d  S  (  <  ( BRP1  +  B RP2) / 2 • 0  >  *  10 • 0 ) 

2  ***  CALCULATE  TH.  cLcCTRICAL  FRcOUENCY  IN  HERTZ 
FRi*  (POL  ES72.) *  (RPH/E.3.3 

2  ***  cAlCJLAT:  THl  electrical  frequency  IN  radians  pep  second 

WF=2. 0*PI*FRE 

2  **•  CalCJLAT.  THPfi  TIMES  TH;.  SKIN  DEPTH  OF  TH.  ENVIRONMENTAL  SHicLJ  IN  IN. 
JELTA=  T. 1»3URT (2. 3 / ( WF»U 2* SI GMA > ) / .  025c 


aoe*Jk-jj?  ,, 
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4  4  4 

♦  4  ■* 


C 

'  *«4 


C  **« 
C  *»♦ 


2  <  * « 
0  *  *  • 
0  *  *  * 
z  *  *  * 


LA.-JL'-T.  U  (JtAL  lUti_".  OK  ACTIVE  CONDUCTOR^  10  TH_  STATjp 
.=(-:♦  j»j)/iph4s: 

;A._Cjl*T-  TH  “r.-OTiVE  .*  C  T I  V  L  STRAIGHT  L  2  N  o  T  H  OF  T  h  l  S  T  A  T  0 
L-PPHER/  (  7. S; -;  o'*  3F  »<;<*«?*•  J*J*:T<G*D) 

CAlCULATl  In.  ACTUAL  PHYSICAL  STRAIGHT  L.NjTm  -ihICH  IS  SMALL;  R  THA..  L  DU  £ 
TO  T-t  .-IF  P>;  HERAT. D  1:1  Trt  lND  TURN  R.GIOIS  OF  T  H_  AIr  GAP  STATOR  WINDING 
lL- L/l .  2  )  1 

CA.CJLAT.  Tn.  TOTAL  .ENgTH  OF  A  RACETRACK  FIclO  CulL  IN  INCHES 
L  C  =  L  L  ♦  -’  •  •  (  J  *  T  ) 

CAlCUlU.  TH.  TOC-  TuD.  1nnE°  rAjIIS  If,  INCHES 
R  i  =  l  {R- A»CuS  (  AMS )  >/-’.) 

CAlCJLA  r.  TH'  TORE  TUI-  UUT:&  •vA  j  1US  IT.  I(OH-_S 

-<  !S  =  R3  +  DlLTAP 

V  IS  TuE  )UT.R  PADIUj  of  the  s  r  A  T  0  3  WINDING 
V-’^-3 

PEG  IS  TnC  OUT-  :  RADIOS  of  TH,  TO.ROU;  TU3E 
R23  =  TRfU.LTAT 

vESO  I s  T Hi  GJTlR  t A  j i us  of  the  :l  .cthomaguetio  shield 
rEjO=r:s+o.ltas 

•!l  IS  Tn-  INNER  RADIUS  OF  T He  TORS  S«- Al 

ri-pis j+g 

<0  IS  TH:  OUTEr  RADIUS  Or  THE  ST  AT  OP  WINDING  (ALSO  EQUAL  TO  V) 

RD  =  R-3  )  *■  o  +•  0  A 

ChLC'JL  A  Tl  TH:.  WIDTH  OF  A  STATOR  COIL  IN  INCH. 3 
U-tO/1,  -»lt) *SORT <1.-C0S(  2.*ANg) ) 

-  IS  TU£  DISTANCE  WHICH  THl  AND  TURN  OF  A  STATOR  COIL  EXTENDS  SEYOND 
TH"  STRAIGHT  SECTION 

.:m/(2*tan(x>) 

CALCULATE  tHF  TOTAL  OV-RAlL  length  of  a  stator  coil 
LPRIMl  =  LL  +  2.*=. 

gALCULATE  THE  COMPONENT  VOLUMES  IN  CU3IC  INCHES  ****♦♦».** 

VOL ( 1 )  IS  THE  VOLUME  OF  THE  STRAIGHT  CYLINDRICAL  PORTION  Or  THE 
.HVIRONMENTAl  SHIELD 


/Ol (1) =PI*L PRIME* (2.*02LTA»RSfD2LT  A**?) 

V0l(2>  IS  TH;  SUM  OF  THl  VOLUME  OF  The  TWO  CUNICmL  portions  of  the 
ln v iron mental  shielc  ^ 


v> 


% 


VO'-  (2)  -  2«*PI*34J-LTA*  (  2.*RSvDELTA-3> 

VOL  (  3 )  IS  THE  VOLUME  Or  THE  TWO  -END  SELLS 
VOL (5) -2. *PI* ( V**2-RRS**2) *OlLTA- 
VO.U)  IS  THt  VOL'JMc  OF  THL  FOUR  FIrL  D  COIlS 
VOl U) =<2.*LL»A*TfPI*A*(2.*B*T*T**2) ) »POlES 
VO- 15)  IS  THl  VOLUME  OF'  ThE  ANTI-DRIVE  END  STRUCTURE 
VO.  (>)  -  0  I  *RRS*  *  2*  (LC/S.I 

VOL(o)  IS  TH;  VOLUME  OF  THE  uRIVE  ENO  STRUCTUR. 

VJl (o)=PI*RRS»*2* ILC/d) 

VOL (7)  IS  THE  VOLUME  OF  ThE  TORDJE  TUSE 
V0L(7>  -F>I»  (2.*0£LTAT*RRtDELTAT»*2>*  (LC  +  S) 

VOL ( 3 )  IS  THE  VOLUME  OF  THE  SORE  TUSE 
VOL  (i>  =PI*LCM2.  *DELTAD*RS+0lLTA3**2> 

VOL(h)  IS  The  VOLUME  OF  THl  fielo  winding  support  structure 
VOL  (  0)  =  (PI’LC*  (RR*»2-R.RS*»2)  )-  VOL  (A) 

VOL  (10)  IS  THE’  VOLUME  OF  THE  ELECT  iOMAGNSTIC  ShIlLO 

VOL  (10)  =  (PI*LC)  *  (RE30**2-P..S**2>  *■  API*  OF L T A 3 > *  ( R£SO**2-RRS»*2) 

VOL (11)  IS  THE  VOLUME  OF  THE  COPPER  STATOR  CONDUCTORS 

VOL (1 1 ) - (2* AC)  * (LL  +W/3IN (X ) ) 


x\%% 


.  b:% 

%y* 


% 


VOl  (12)  IS  THl  VOLUME  JF  THE  STATOR  INSULATION 
VOL (12) = (PI*LPRIM2>* ( RO* » 2-KI* *2 ) -VOL (11) 


O  C>  O 
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C  ***  VOL  (15)  la  THE  VOLUME  DF  THfc  BORE  SEAL 
VOL (13) =PI»LFRIME*  (2.*  OUSEL « R  + 09SEL *  * 2 > 

0  **»  Qt.c'IM1:  THE  COMPONENT  MASS  DENSITIES  in  POUNDS  PlR  cubic  inch  *** 
p< 1)  =  .  1 
P (2 ) -•  1 
P  (  3 )  =  • 2  8 
P  (  w )  = .  2  G  5 
P ( 5 ) = • 23 
P  ( 6 )  -  •  2 
P(7>=.23 
P( 3>=. 23 
P<  )>=.  1 
P(10)=.l 
P(ll>=. 32 
P(12)  =.  1 
P  ( 1 3)  =  •  10 
WT - C  •  3 

C  ***  CALCULATE  THl  weight  of  the  rotor  in  pounds 
ROTWT-j.  0 
DO  43  J  =  4  » 1  0 

43  kOTWT=ROTWT+P(J)*VOL(J) 

C  ♦»»  CALCULATE  THl  TOTAL  ENVELOPE  VOLUME  OF  THE  SUPERCONDUCTING  GENERATOR  ***** 
C  *»*  THE  units  are  in  cubic  feet  ******* 

VOLUM£=(  LPRIME*PI*(RS*OELTA) **2  *  2.*S*PI»(  2  .  *RS**2 «■  S  **  2-2.  * S*  RS 
4 »RS* ( RS-S)  )  )/1728.0 

C  »*»  LALCULATc.  THE  TRANSPORT  CURRENT  LOSScS  (OH-IIC  HlATING)  IN  KILOWATTS  *** 

C  ***  CONVERT  THE  COPPER  VOLUME  IN  CUBIC  INCHES  INTO  CUBIC  METERS 
VOLCU=VOL <111*1. 6387E- 05 

C  **»  CONVERT  THE  CONDUCTOR  CURRlNT  DENSITY  IN  AMPS/SQ. IN.  INTO  AMPS/SQ. METER 
A JO  =  JC* l . 55E  +  0  3 

E  **♦  GA..CULATS  THE  RESISTIVITY  OF  COPPERfRHO  IN  OHM-METERS)  AT  TMPC  DEGREES  C  * 

T M3C  =  1 5 0  •  0 

RHO=l. 72E-C3* (23m.5*THPC) /254.B 
C  ***  PI RL  IS  I  HE  TOTAL  I  SQUARED  R  LOSS  IN  KILOWATTS 
PIRL=AJC*AJC*RHO*VOLCU/lCOO.uO 

C  *»*  CAlCULATl  THE  EOOY  CURRENT  LOSSES  IN  THE  STATOR  BARS  IN  KILOWATTS  »♦* 

C  »**  CALCULATE  THfc  PEAK  TANGENTIAL  MAGNETIC  FLUX  UENalTY  AT  BOTH  THE  INNER 
C  AND  OUTER  RADII  OF  THE  STATOR  WINDING  AT  AN  ANGLE  PI/4  (IE  HALF  WAY 

C  id W_£N  POLES)  FROM  THE  COIL  BORE  (SEE  FIGURE  30  IN  THIS  APPENDIX) 

PHD- j . 25  j 

CAlL  MCLD(RM1,PMU,9R,3TM,X0 , Y 3 , OX  ,  OY , A I ) 
call  MFL J<  RM2.PMU, BRR, 3FF,XC,YD,DX,DY,AI) 

3  ***  CALCULATE  THE  AVERAGE  OF  THE  TWO  P-.AK  VALUES  OF  RADIAL  MAGNETIC  FLUX 

C  ***  DENSITY  AT  THE  INNER  AND  OUTER  STATOR  WINDING  RADII 

oRPM= (  )R1kBR2)/2. D 

C  ♦**  CALCULATE  THE  AVERAGE  OF  THE  TWO  PEAK  VALUES  OF  AZIMUTHAL  MAGNETIC  FLUX 

0  *«*  JETSiTY  AT  THE  INNER  AND  OUTlR  STATOR  WINDING  RADII 

ATPM- ( dTM»PFF)/2, 3 

C  ***  CALCULATE  THE  EJJY  CURRENT  LOSSES  IN  WA TTS/CU3 IC  METER 

PE/=(WF**2/ <2-t.*RHO) )* ( (H0**4-TI**-*) / ( MC ** 2- T I ** 2 ) ) * ( 9 TPM**2+ 3RPM* 

4*2) 

*»♦  HI  )UCE  THE  c.FFlCT IVE  VOLUME  OF  STATOR  COPP:R  TO  ACCOUNT  FOR  THl  FACT  THAT 
*»*  HOST  OF  THl  LJJY  CURRENT  LOSSES  OCCUR  IN  THE  STRAIGHT  SECTION  OF  THE 
***  ST  A  T  JR  WINDING 

VOICE =VOLCU*LL/(LL*W/SIN (X)  ) 

0  **♦  LAlCULATl  The  total  EOjY  current  loss  in  KILOWATTS 
P-_  1D=P-V*VOLCE/1g  JC.uD 
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***  oAlCJLAT-E  THl  LOSSESIEDDY  CURRENT)  in  the  image  s'hield  in  kilowatts 
♦**  del  is  the  skin  depth  in  the  aluminum  fnvironmental  shield 

0EL=SURT (2. / (WF*UO»SIpMA)) 

i<Sv(  r9  THE  RADIUS  OF  THl  ENVIRONMENTAL  SHIELU  IN  METERS 
RS  i=RS* • 0254 

»♦*  LSM  IS  THE  LENGTH  OF  THE  CYLINDRICAL  PORTION  of  the  ENVIRONMENTAL  SFiIElO 
♦»*  IN  METERS 

LSM=LPPIMZ*.C25N 

**»  CALCULATE  THI  PEAK  VAlJE  OF  TH_  A7 IMUTHAL  COMPONENT  OF  MAGNETIC  FuUX 
**♦  DENSITY  AT  RADIUS  RSM.  NOTE  THAT  THE  PEAK  OF  AZIMUTHAL  FLUX  DENSITY 
*»*  OCCURS  MID  WAY  NET  WEEN  POLl  CENTERS,  THAT  IS,  AT  AN  ANGLE  OF  PI/A  FROM 
***  THE  COIL  30RE  (  SEE  FIGURE  30  IN  THIS  APPENOIX  ) 

PMO=D. 250 

CALL  MFl  MRSM,PMU,BP, 1T,XO,YC,OX,OY,AI) 

***  PSHLO  IS  THE  EOOY  CURRENT  LOSS  IN  THE  ENVIRONMENTAL  SHIELD  IN  KILOWATTS 
PSNlD=  <PI*RSm»lSM*GT*3T> / (aIGMA*UEL*U0*U0*lO00.C3) 

calculate  the  generator  efficiency  (percent) 

EFF= (POWER/ (POWER* 100 u.C» (PIRLyPEDD+PSHlJ) )) *1jC. 00 
*•*  CONVERT  THE  CONDUCTOR  HEIGHT  FROM  METERS  TD  INCHtS 
HOI=HO»1CO. 0/2. 5a 

***  CONVERT  THE  CONDUCTOR  WALL  THICKNESS  FPCM  METERS  TO  INCHES 
TwI=TW*l J0, C/C.Sh 

***  CALCULATE  Th.  POWER  IN  MEGAWATTS 
PM=POWER»1.E-06 

*•*  CONVERT  T He  END  TJRN  EXTENSION  ANGLE  FROM  RAOIANS  TO  DEGREES 
XA  )= ( X/PI ) » 10b. 0 
PRINT  d0 

PRINT  73,  PM,VJC,RPM,FRE, JC , V T , PF , F9 AR , KW 
PRINT  91 

JRI NT  7 1  ,CS, ANS,  Z, IPHA SE ,  SLO T S , SW I UTH , DO AR, OCD N , GPH 
PRINT  92 

POINT  72,  DElTAU, JLLTAT, DELTAS, 0E L f AF , S ,G , D , A, XAN 
PRINT  33 

PRINT  73,  RR, D ,RS , LL ,L C, lPRI Ml , E , O JSE L , T 
PRINT  i, 

PRINT  '7,  JA, D ARM, H0I,TWI,W, DELTA, FARM, AC 
PRINT  35 

PRINT  7<*,  'JR1,HR2,GKG,  )RM,rlTM,  3FF.uRPM,f3TPM,BT 
DO  53  <=1,13 
nT  D-P ( < ) "VOL (K) 

PRINT  j? t  K  ,  PTC 
5"  WT=WT+NTC 

***  CAuCJLAT_  THE  SPECIFIC  W.IGHT  IN  POUNDS  PER  KILOWATT 
jwt=io; d.*wt/powzr 
PRINT  36 

PRINT  ’5,  SWT, WT.POTWT, VOLUME, Q, ACM 
PRINT  3? 

PRINT  76,  EFF,PIRl,PEDU,PSHLu 
331  GO  TJ  1 5  j 
IS!  PRINT  13* 

153  CONTINUE 
b~  FORMAT  (UJ,F12.2) 

7:  FORMAT (9F12. 2, /> 


7 1  FORMAT (3F12. 2, /» 

72  FO<MAT ( 5F12. 2,/) 
y3  FORMAT < 9F12. 2, /> 
7 h  FORMAT ( 3F12.E,/) 


THIS  WCV  iKf'T  IS  mF 


,  .  rP7i.ETISEifiI«E. 


r-.R 

JHGNIPICANT!  KUK3ER  0d?  FACSS 


■  ■Hf’ilOJUiCiL  lESGLfcli*  °i 
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FORMAT  (F15*4,5F12.2,/) 

7o  FQ >MAT (F15. A,  3F12.2) 

77  FORMAT (8F12. 3, /> 

8  9  FORMAT  (  lnl ,  hX^HPIMH)  ,  7X  ,  8H0C  VOLTS.oX  ,  3HRPM  ,  7  X  ,  “.HFREQ  ,dX,2HJC,12X 

A,  2  Mi/ 1  ,  1 JX.2MPF, 3X,4HF3AR,  luX, 2HKW) 

9  1  FORMAT  !«X,2HCS,UX  ,2rtNS,  10  X,  1HZ,9X,  3HIPH,  7X  ,  5HSLOT3,7X  ,  6HSRITTH,  9X 

i,-*HGjAR,7X*4HOCON,10X,  3HPPH) 

3  2  Fo  RMAT  <7X,-*HDELU,  8X,4M0ELTf  9X  ,  4  H'J  t  L  S  , 8 X , 4H TEL E , 1 C X , 1H$ , 13  X,  3rl&AP  , 1 
ilX,lr!J, 11X,1HA,9X,7HE JO  AN&) 

9  3  FORMmT l 3X , 2HkR,UX,  1H3, 10X  ,2HRS, 1 C X , 2HL L , 9X , 2HLC , 9X , oHLPR I  ME,  9X , 1H 
i_,  3X,5HDdSSL,9X,lHT> 

3  4  PO  RMAT  l 8X, 2HCA ,8X,4H0ARM  ,  1 IX  ,  1HH  ,  1 3  X  ,  2wT  W ,  1 0  X  ,  lh  A ,  MX  ,  5H0.-.LT  A ,  7  X  ,  4M 
LFA RM, 8X  , aHACOND) 

9  3  F JRMAT ! 9x, 3H  TR1, MX, 3H3R2,  9X  ,  3rl3Kr> ,  MX  ,  3H3RM  ,  9X,  jFmGTM,  3X  ,  3rf3FF ,  8  X  ,  4rt 
i  ORPM ,  7  X  ,  ■* H 3  T PM ,  1 C  X  ,  2 F!  3 T ) 

<So  FORMAT  l  /  ,  2X  ,  1 3rtSP£C.  WT .  L  o/ KW  ,2X,1QMTOT.WT.lBS,3X,  3HROT  OR  (  L3)  ,  2X, 1 j 
iMV'JLtrJ.FT)  ,2Xi1JHSPEC.  lOAO) 

97  FORMAT (BX,1 3HEFFICIEM3T, 3 X , 1 fl HOHM 13  LOSS , 3X , yHiD^T  L JS 0 , 2 X , 11 HSM IE 
LL3  LOSS) 

1 F 1  FORMAT!//, 13h  NO  ROOM  FQR  COPPER* / ) 

ST  J3 
ENT 
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SO  JROUT  Xfj;  1FLu<  C-IO,P  TU,  ORHO.dPHI  ,Xt. ,  Y*  ,OX, OY  ,  61) 

JO  J  1...  PRECISION  -t,d,AI,  X  0  ,X1 ,  YS  ,DX,  0Y,PI,P2,AM,'K,UF  ,XXT  ,  YYC 
00  J  ml  '  Ok;GIjI  J'l  XX,  YY  , RH  0 , Z  ,  PH  I ,  X,Y  ,  F ,  DKHO  i  OPrt  I ,  QlAG 
Jl*  MS  TON  AC.2,12,J),j<12,12,>i) 

^  ♦  ♦♦*♦**♦♦♦#•*♦**♦ 

C  •»*  HIS  SIUTOUTIM;  CuLC’JLITf.O  T Mt  MA&M-.TIC  FLOX  OLNSITY  IN  FRI-  SPxO£ 

0  ***  F0J.  A  FOUR  POL'  KOTOR  3TJUCTUR-  JSINS  T  HI  UOT-SAMART  LAN. 

0  ***  A  12X1  ?  ARRAY  OF  MISTS  Rt_i-RI  S£  NT  S  TrtCH  COIL  SI  O'-  WriICrt  13  ASSUNiO  TO  1-- 

o  *»*  rectangular  in  cross  section. 

0  •••  Tri  RECTANGULAR  jTCSS  SICTION  IS  (12*P X)/.025j  INCli.3  NIjE  JY 
3  ***  (12*0Y  >  /  .  C2S-.  INCH  S  HIGH.  wh£R'  HIJc  AND  HIGH  RlFlR  TO  ToS  X  AN  J  Y 
C  **♦  C01RTMAT:  01-UCTICMS  R£SP; CTIMELY.  (SSL  F IGuR-  30  IN  THIj  APPENDIX) 

^  ♦*♦,»»****»•*♦,♦ 

P I  =  3 .  I<l5i2u53-5391 
JH  T  =  DN  J*PI 
°2  =  P I / _  . 

.X  =  |.,  i-C7  *PI 
CR=A  1»  11/  (2,  »“I) 
x*=-x: 

K  =  1 

<  =  X  C 

Y-Y  j 

3'  CO  1  J=1 , 12 
CJ  2  1=1,1? 

- < : , J,<) =0SQKT (X**2mY*»2) 

7=  Y/X 

111,1,0  =JATAN(Z) 

2  X  =  X  M  ix 

x=  < .. 

1  Y-Y ♦ J  Y 
<=2 
X  =  X1 
Y=YT 

JO  '  0  =  1,12 
00  L  T  =  1 , 1 2 

1(1,  J,<)  -  )SNPT  (X-»»2  +  Y**2) 

Z=  Y/X 

-  (  [  ,  J,  <  >  =  0A  T  AN ( Z ) 

-  X=X-JX 
X  =  <  1 

0  Y  =  Y  +JY 

<  =  3 

00  5  J  =  1 , 1 2 
CO  5  1-1,12 
A  ( I  ,J,  X)  =  A  ( I  ,  J  ,  1 ) 

3  3(1, J, 0=0(1, J,l)+P2 

<=, 

00  t  J=l,12 
00  6  1=1,12 
A(I,J,<) =A ( I ,  J , 2) 
o  y (I , J, K) =0(1 , J,2) +  P2 
K-i 

00  7  0=1,12 
00  7  1=1,12 
a ( i ,  j,o  =Ad,  j,  3) 

7  3(I,J,<I =3(1, J, i) +P2 

K  =  ■, 
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10  6  J=l,12 

00  <3  1=1,1? 

A(I,J,<)=A(I,J,4> 

3  11 I , J,K>  =  3(1  ,  J,4> *P2 
*=7 

01  9  J-  1,  12 
00  9  1=1,12 
A(I,J,0=A(I,J,5> 

0  0( I ,  J, K) =  9( I , J,3) +°d 
K=3 

00  i:  J  =  1 » 1 2 
00  11  1=1,12 
A(I  ,  J,<) =A  (I ,  J,  a) 

10  9(1  ,J,K)  =9<I,  J,6n-P2 

mno=a.  oo 

9PrlI  =  0  .00 
jF  =  -  3< 

K=1 

00  11  3=1,12 

DO  11  1=1,12 

F  =  A  (I,  3,K>  **2+9H0**2-2.  ♦  A  (I ,  3,9)  ♦  P.,90  *UC0i  <3  ( I ,  J  ,  K )  -Pri  X  ) 
■3990=99H0t-bF*A  (I  ,  3,9)  ♦OSINIB  (I  ,  J  ,  K 1  -PHI)  /F 

11  OP  31=  3DHI<-9F«  (PHO-A ( I,  3, 9 >  *DC03 ( 3 ( I , J , 9 ) -P9 1 ) > / F 
9F=+39 

9=2 

00  12  3=1,12 
00  12  1=1,12 

F  =  A (I,  J,K) **2«-RHU*»2-2.*A ( I ,  J  ,  9)  * RrtO*  3C03  (  J ( I , J  ,K) -PHI ) 
9RH0=9PH03-9FM(I,  J,K)  *0SIN(B(I,  J,K)-PHI)/F 

12  JP9I  =  0PHIH0F*  (PHO-A  (Ii  J,K)*DCOSU(I,  J,K)-PHI)  )/F 
9=3 

00  13  J=  1 » 1 2 
00  13  1=1,12 

F=4 (I, 3,9)  »*2+PHO**2-2.*»  <1 , 3,9) »PH0*DC03 <  3 < I ,  J , 9 ) -PH  I ) 
3PH0=09H0»9F*A (I , J , K > *0S IN < 6 ( I , J , K ) -PHI ) /F 

13  OPHI=3jHIh'3F*(PHO-A(I,J,K) *DC0S(9(I,J,K>-PHI>)/F 
3F  =  -  3K 

9=4 

00  1-4  J  =  1 , 1 2 
DO  14  1=1,12 

F=Ul,3,9)**24RH0»*2-2.»A  ( I ,  J ,  9)  *  PHO*  JCOS  <  3  ( I ,  J,  <)  -PHI ) 
3990=9  PH0*3F*A(I,  J,  <>  ♦  3S1M  ( *3  <  I  ,J,K)  -PHI)/F 

14  JPH  I=I3PH  I +  9F*  (PHO-A  (I,  J,K)  ♦DC0  3<  3(1  ,  J,K)  -PHI)  )  /F 
9=  3 

00  15  J= 1 , 1 2 
00  15  1=1,12 

F  =  A(I,J,<)  **24RH0»*2-2.»  A  (I,  3,9)  *PHO*OCOi(  3  ( I ,  J  ,  9) -PH  I ) 
3  99  0=  999 J*  9F  *  A ( I , J , 9) * OSI N ( 3 ( I , J , K) -PHI ) /F 

13  JP9I  =  3PHI«-3F*  (PHO-A(I,  J,K)  *QC0S(3(I,J»K)-P9I))  /  F 
bF  =  ♦  39 

9=5 

00  16  3=1,12 

00  lo  1=1,12 

F  =  A (I, J,K) **2+PH0**2-2.»A (I, J,9) *RHO*0COS< 3 ( I , J , K ) -PH  I ) 
rpH0=0  PH090F*A(I,J,K)*0SIt((8(I,J,K)-PHI)/F 

14  0PHI=9PH I F3F* (PHO-A (I, J,<) ♦OCOS( 3( I, J,9)-PHI) ) /F 
9=7 


.......  -  —  f*RACTIC«W«*Ji 

•  *  .  ■  ?  r  g  -  ~  v 

X 

it .  1  ■  ■  ;  rfl  Mia  f 

>XUr  :  "  -V  . 

•4<i’  4'  *.(&  UDGIidJjXo) 
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30  17  3=1*12 
00  17  1=1,12 

F  =  a  <  I ,  J,<)*»2«-RH0*»2-2.»&(I,J,O»Rrl0«QC0S<  HI,  J,K)-PHI) 
RHO=  V7HO*3F*A (I , J,K> ♦OSIN (0 (I , J,K> -PHI) /F 
17  dP-tl^^HI  +  OF*  (P.HO-A(I,  J,  O *DC0S(9(I, J,K>-PHI) > /F 
3F=-3< 

K=9 

)0  lo  J=  1, 12 
JO  19  1=1,12 

F=  \  (I,  J,<)  *  *  2  +  RHO*  *  2-2 .  *  A  (1,3,0  *KrlO*OCOS  (  M  I,  J,K)  -PHI ) 
3-H0=  3  <H07  JF*4(I,JfO*JS  I  N ( B ( I , J , K> -PHI) /F 
13  3°9I=9^HIt3f* (RrtO-MI ,  J,  K. ) *DCOS O ( I , J , K ) -PH  I > ) /F 
PtTUP.N 
INO 


SAMPLE  OUTPUT 

SUPERCONDUCTING  GENERATOR  BASE-LINE  DESIGN 
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APPENDIX  C 

POWER  SUPPLY  COMPONENT  WEIGHT  ALGORITHMS 

The  following  weight  algorithms  were  used  to  calculate  component 
weights  for  the  system  point  designs  described  in  Table  6. 

FUEL 

WpuEL  =  0.3724  (SPC)  PjURB(M)  AtQn  (lbs) 

SPC  =  specific  propellant  consumption,  1 bs / ( HP-HR ) 

PtuRb  =  turbine  shaft  power,  megawatts 

AtQn  =  total  power  generation  time,  seconds  (assumed  to  be  120  seconds 
for  the  point  designs) 

TURBINE  (Reference  2) 

^TURBINE  =  0,3455  °T  *  (1bs^ 

where 

Dt  =  3.896  x  lO^/RPM-p  (Turbine  tip  diameter,  inches) 

RPMT  =  Turbine  rotational  speed 

GEAR  BOX  (Reference  11),  Double  reduction,  three-branch  gear  box 
WQB  =  (88.236/RPM-j.)  (HP)  (SUM),  (lbs) 

where: 

HP  =  mechanical  shaft  power  of  the  turbine  in  horsepower 
SUM  -  }  0  ♦  mg)  ♦  mg  (2  ♦  mg  *  »g/H0)  «  ^  M0  0  +  tyl^) 

M„  =  RPMT/RPMr  =  overall  gear  ratio 

0  I  b 

RPMg  =  rotational  speed  of  the  generator  rotor 
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m^  is  the  optimum  internal  gear  ratio  and  is  found  from  the  following 
expression : 


2m2  [m  (M  +  1 )  +  M  ]  =  i  M  (M2  +  1 ) 
g  L  g  o  '  o  3  o  o 


A  linear  approximation  to  the  above  expression  is  given  by  (R.  L.  Binsley, 
Rocketdyne,  personal  communication): 


M  =  0.185  M  +  .333 
g  o 


For  each  of  the  point  designs  in  Table  6  involving  a  gear  box,  Table  9 
lists  the  pertinent  parameters  for  calculating  the  gear  box  weights. 

The  weights  were  calculated  for  both  the  exact  and  approximate  values  of 
m^;  however,  the  exact  values  of  m^  were  used  in  calculating  the  gear 
box  weights  in  Table  6. 


TABLE  9 

GEAR  BOX  WEIGHT  CALCULATIONS 


POINT 

DESIGN 

rpmt 

rpmg 

M 

0 

mg 

(approx) 

(exlct) 

WGB(lbs) 

(approx) 

WcB(lbs) 

(exact) 

2  ( PMG) 

14,900 

18,000 

0.82778 

0.48614 

0.38880 

204 

182 

3  (SCG) 

14,900 

6,000 

2.48333 

0.79272 

0.76030 

487 

484 

5  (PMG) 

10,500 

18,000 

0.58333 

0.44092 

0.340675 

602 

507 

7  (SCG) 

10,500 

6,000 

1  .7500 

0.65675 

0.592207 

1183 

1156 

HIGH  POWER  TRANSFORMER  (Tl),  (References  4,  13) 


W 


Tl 


0.931  +  0.069 


fe)"]  ■  [•■ 


242  +  0.758  f 


-0.926 


X  P 


1(KW) 


(lbs) 
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VQ  =  DC  output  voltage  required  at  the  load  (KV) 
f  =  alternator  frequency  (KHz) 

P1(M)  =  transformer  output  power  (MW) 

P1(KW)  =  transformer  output  power  (KW) 

LOW  POWER  TRANSFORMER  (T2).  (References  4,  13) 


0.612  +  0.388 
0.71  1 


-0.985 


X  f-°-767 


X  P2(KW) 


(lbs) 


P2(M)  =  outPut  P°wer  of  low  power  transformer  (MW) 
P2(KW)  =  outPut  Power  of  low  power  transformer  (KW) 

THREE  PHASE  RECTIFIER  (References  4,  13) 


Wrect  =  °-0073 


0.945  +  0.055  ttttt 
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X  P(KW)  <'bs) 


P(KW)  is  tbe  ’nPut  power  to  the  rectifier  in  kilowatts  and  is  P](kw) 

for  the  high  power  branch  of  the  circuit  or  P2(kw)  for  the  low 
power  branch. 


NOTE:  VQ  must  be  chosen  correctly  to  reflect  the  correct  weights  for 
the  high  and  low  power  branches  since  VQ  is  different  in  the  two 
cases . 


LC  FILTER  (References  4,  13) 


"LC 


1.81  x  10 


-3 


X  P( KW)  »  (lbs) 


89 


AFAPL-TR-79-2073 


where  p(KWj  is  defined  under  THREE  PHASE  RECTIFIER.  The  following 
assumptions  were  made  concerning  the  LC  filter: 

1%  peak-to-peak  ripple  allowed  in  the  load  voltage; 
inductor  energy  storage  density  =  25  joules/lb; 
capacitor  energy  storage  density  =  100  joules/lb; 
conduction  overlap  angle  =  30  degrees; 
diode  rectification  (i.e.,  no  firing  angle  control). 

PMG  COOLING  SUBSYSTEM  (Reference  2) 


W 


c 


130  P 


0.375 

G(M) 


(100  -  nG) 


X  At. 


on 


(lbs) 


PG(M)  =  9enerator  output  power  in  megawatts  times  the  number  of 

generator  units  required  to  supply  the  required  load  power 

nG  =  generator  efficiency  in  percent 


SUPERCONDUCTING  GENERATOR  COOLING  SUBSYSTEMS 
Stator  Cooling  (Pressurized  Water) 

000  -  n  ) 

W„  =  x  - —  X  W„ 

c  G(M)  nG  c 

Wc  =1.0  lb/(MW-SEC) 


Cryogenic  Cooling  (Liquid  Helium) 

Wcryo  =  (l.°9)  n.l°)  (L)  (MT)  (lbs) 

1.09  =  0.80  +  0.29  accounts  for  the  weight  of  the  dewar  system  and 
the  stored  liquid  helium,  respectively. 

1.10  =  factor  to  account  for  losses 

L  =  rotor  helium  flow  rate  required  in  liters  per  hour.  For  the  10  MW 

machines,  L  is  assumed  to  be  20  liters/hr  and  for  the  20  MW 

machines,  L  is  assumed  to  be  25  liters/hr. 

MT  =  mission  time  (duration)  in  hours,  assumed  to  be  10  hours. 
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